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Abstract 

These lectures review the progress made in our present understanding of B decays. The 
emphasis here is on apphcations of QCD to B decays and the attendant perturbative and 
non-perturbative uncertainties, which hmit present theoretical precision in some cases 
but the overall picture that emerges is consistent with the standard model (SM). This is 
illustrated by quantitatively analyzing some of the key measurements in B physics. These 
lectures are divided in five parts. In the first part, the Kobayashi-Maskawa generalization 
of the Cabibbo-GIM matrix for quark flavour mixing is discussed. In the second part, 
the bulk properties of B decays, such as the inclusive decay rates, semileptonic branching 
ratios, 5-hadron lifetimes, and the so-called charm counting in B decays are taken up. 
The third part is devoted to theoretical studies of rare B decays, in particular the elec- 
tromagnetic penguins involving the decays B —>■ K* + 7 and B —>■ Xg + j. The photon 
energy spectrum and the branching ratios in the SM are discussed and compared with 
data, enabling a determination of the CKM matrix element \Vts\, the 6-quark mass, and 
the kinetic energy of the 6-quark in the B meson. The CKM-suppressed inclusive decay 
B — > Xd + 7, and the exclusive decays B (p, a;) -|- 7, are discussed in the SM using 
QCD sum rules for the latter. The importance of these decays in determining the CKM 
parameters is emphasized. This part ends with a discussion of the decay rates and dis- 
tributions in _B — > Xgi'^i", including the long-distance effects, and estimates of a large 
number of other rare B decays, including B — > Xuu, and the two-body decays B'^ — > 
and B^ — *■ 77. The fourth part is devoted to reviewing the present estimates of the CKM 
matrix elements from B decays and B^ - B'^ mixings, which determine five of the nine 
elements in this matrix. This is combined with our present knowledge of the other four 
CKM matrix elements and a quantitative test of the unitarity of the CKM matrix is 
presented. This information is then combined with the constraints from the CP- violation 
parameter |e| in order to provide a profile of the CKM unitarity triangle and CP- violating 
asymmetries in B decays. These aspects are discussed in the fifth part of these lectures. 



* Lectures given at the XX International Nathiagali Summer College on Physics and Contemporary 
Needs, Bhurban, Pakistan, June 24 - July 13, 1995; to appear in the Proceedings (Nova Science Publishers, 
New York), Editors: Riazuddin, K.A. Shoaib et al. 
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1 Weak-coupling Lagrangian for the b Quark and the 
CKM Matrix 



The Lagrangian density of the standard model for electroweak interactions (henceforth 
called SM) can be symbolically written as JIl: 

C = C{f, W, B) + £(/, $) + CiW, B, $) - \/($) (1) 

where the symbols /, W, B and $ represent fermions (leptons and quarks), SU (2) ^ gauge 
bosons, W^, U{1)y gauge boson, and the Higgs doublet field, respectively. The SM 
particle content together with the (weak) isospin properties of the basic quanta are given 
below, where we have assumed that there are three families of leptons and quarks. 



Fermions 



Leptons: 






Gauge bosons 



Scalars 



- /AO 



(2) 



The various terms in the SM Lagrangian can be written by demanding SU{2)l U{1) 
gauge invariance, lepton-quark universality, and family-independence of the electroweak 
interactions. Since the SM Lagrangian is given in any standard textbook on electroweak 
interactions we shall not reproduce it here. The main interest in B physics lies in the 
study of the first two terms in (1) involving fermions, gauge bosons and the Higgs fields, 
in particular flavour-changing transitions. In order to appreciate how flavour- changing 
transitions emerge in the SM, it is worth while to write the first two terms in (1) explicitly. 
The interaction between the fermions and gauge bosons has the form: 



(3) 



2 



where j is the family index, 



^i=(^g^j , Ir = eR , ql= (^^^ ,q]i = UR, ql^dR,... 
and the two covariant derivatives are defined as 

Di^^d^- ig2 (w^ ■ - igi^Bf,, 

D'^^d,-ig,^B„ (4) 

where gi and g2 are, respectively, the U{1)y and SU{2)l coupling constants, a" (a = 
1, 2, 3) are the (weak) isospin Pauli matrices, and the Gell-Mann - Nishijima formula 
Q = Is + Y defines the (weak) hypercharge of the quarks and leptons. The interaction 
term £(/, 0) involving the fermions and the Higgs fields has the Yukawa form 



>c(/,$) = 

3=1 

+ E { {K),, d^^R + ih.hkqi^'d'k], (5) 



where 



j,k=i 



with both $ and transforming as a (weak) isospin doublet. The hypercharge of the 
Higgs fields can be written by inspection. The Yukawa coupling constants {hi)j, {hq)jk 
and {h'gjjk are arbitrary complex numbers and each term in (5) is independently SU{2)l® 
U{1)y invariant due to the fact that the SU(2)l acts only on 1^ and qL and on the Higgs 
doublet $ and whose products are SU (2) scalars. 

After spontaneous symmetry breaking SU{2)l ® U{^)y — ^ U{^)em-, which preserves 
the symmetry under ^7(1) electromagnetism, the gauge bosons, fermions and the neutral 
scalar field acquire non-zero masses through the Higgs mechanism 

l-(0) = //2|0|2 + A|0r ; //'<0,A>0, (6) 



with y (0)min at |0| = v/ \/2 = y — /x^/2A, where v is the neutral-Higgs vacuum expectation 
value. Making the Higgs transformation (f) — > (f) + v in (5), one finds {(f) is the physical 
Higgs scalar with = 2Av^) : 



7 = 1 ^ ^ 



3 / 1 \ 

E {Kfe)[/44 + Kfc)D 44} 1 + -0 + h. c, (7) 

i,k=l ^ ^ ^ 
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where 



Since in the SM there are no right-handed fields z/)^ {i = 1,2,3), the neutrinos re- 
main massless and the charged lepton mass matrix {mj)^ is diagonal. Hence, in the SM 
there are no family- changing leptonic interactions - an aspect that will soon come un- 
der experimental scrutiny from the ongoing and planned experiments dedicated to the 
neutrino mass and oscillation measurements. In (8), {mjk)^ and {mjk)j-) are the (3 x 3) 
quark mass matrices for the up- and down-type quarks, respectively. In order to write 
the Lagrangian in terms of the quark mass eigenstates one has to diagonalize the mass 
matrices {mjkjjj and {mjk)jj. This can be done with the help of two unitary matrices. It 
is customary to denote them by V^^a,nd V^^^ (likewise for the down type quarks): 

Vl^^muV^'^^ = (mdiag.)^ = Diag. (m„ ,mc ,mt) , (9) 
V^°"°mi:, V^°"°^ = (mdiag.)^ = Diag. (m^ , , m(,) , 

with V^P^V^P = 1, etc. Concentrating on the up-type quarks in (7) one can do the 
following manipulation : 

uluiuUr = ULV^^W^^muVl^'^Wl^^UR (10) 

= Vl'^Ul (mdiag.)^/ {Vr^Ur) , 

which shows that the physical quark fields (mass eigenstates) are: 

^^_(Phys) ^ (^ydo^n^^y_ ^^^^ 

Likewise, u^^^^'' = {V^^urY , (t^^^^^ — (Vr""'^^^) . One can now rewrite the term 
£(/, $) in the SM Lagrangian in terms of (u^^)^^^^ and (d^)^*^^^, obtaining 

C{f,<^r''^-(^l + fj {Em,,g-g, + |:m,/,/,}, (12) 

where it is now understood that ?i = | (u^^^^ + u^^^^^^ = u etc., and we have dropped 
the superscript on the quark fields. The identification or the parameters mi^, mq^ with the 
lepton and quark masses is now evident. In addition, the SM Lagrangian has specified 
the Higgs-fermion (</>//) Yukawa couplings and their Lorentz structure, as well as their 
C, P, and CP properties. All these symmetries arc separately conserved in £(/, 0)^^^ and 
the Higgs-fermion Yukawa couplings are manifestly diagonal in flavour space. It should 
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be emphasized that this is a consequence of the choice made in the SM of a single doublet 
Higgs field since, otherwise, flavour-changing neutral-current (FCNC) transitions in the 
Higgs sector would be allowed in general. 

Finally, one can carry through the transformation (10) in the part of the SM La- 
grangian describing the fermion-gauge boson couplings, C{f,W,B). Written in terms of 
the physical boson (PF^, Z'^, A^) and fermion fields, it is easy to show that the neutral 
current (NC) part of C{f, W, B) is manifestly flavour-diagonal. Thus, all flavour- changing 
transitions in the SM are confined to the the charged current (CC) sector. Denoting the 
quarks and leptons by fiii = 1...6), the neutral current in the SM is given by: 



JH"" = E \ ■ n ^ a (^3L - Q Sin ''Qw) . + eA^Q, 
^ ^ sm t^VK cos b'H/ ^ ^« 



(13) 



where = (1 ~ 75)/2(/3) with J3 = +1/2 for Ui and z/j and —1/2 for di and Qi is 

the electric charge of the fermion /j in units of the electron charge, i.e., Qf. = +1. The 
electroweak mixing angle in J^*^, denoted by 6\y, has its origin in the diagonalization 



of the gauge boson mass matrix, and it has the usual definition cos6'iy = 92/^91 + 9^ 
with the electric charge defined as e = 92 sin 6 w It is easy to check that the neutral 
current interaction induced by the Z exchange violates P and C but conserves CP. The 
electromagnetic interaction conserves, of course, all three C, P, and CP separately. 

The neutral current couplings specified in J^*^ have been measured in e~^e~ annihila- 
tion experiments, in particular at LEP and SLC, and elsewhere. A comparison of these 
couplings with data can be seen in the Particle Data Group (PDG) review of particle 
properties 0. The present situation can be summarized by the statement that the SM 
is consistent with the vast majority of these measurements. Quantitatively, this can be 
expressed in terms of the following values of the basic parameters of the SM, obtained by 
averaging the LEP, SLD, pp and uN data 

mz (GeV) = 91.184 ± 0.0022, 

mw (GeV) = 80.26 ±0.16 

sin^ Owimz) = 0.2315 ± 0.0002 ± 0.0003 , 

mt (GeV) = 180 ± 7^^^ 

as{mz) = 0.123 ± 0.004 ± 0.002 , (14) 

where sin^ 9 w{mz) is defined in the modified minimal subtraction scheme (MS) at the 
scale mz, and the particular definition of sin^ 9]y used in the analysis can be seen in 
Langacker's review in [Q. In addition, N^, = 2.991 ± 0.016 0], where N^, is the number of 
light neutrini. The fitted value of mt from the electroweak data is in excellent agreement 
with the direct measurements of the same at the Fermilab experiments CDF and DO 
mt = 180 ± 12 GeV. However, precision measurements at the Z^ also yield decay rates 
r{Z bb) and (to a smaller extent) r{Z — > cc), which are on their face value at variance 
with those predicted in the SM, taken together with the measurement of the top quark 

tThe numbers are taken from P. Langacker's 1995 updated review available through the World Wide 
Web. 
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mass. The experimental situation at the end of 1995 has been summarized in [0, 



^'"f(Slsb = "■2™ ±0,0017 |ft(SM)= 0,2156], 

^ tedrons) = 01543 ± 0,0074 |fl,(SM) = 0,1724], (15) 

These measurements constitute a pull factor of 3.7 and —2.5 on the SM values of Rb and 
Rc, respectively ^. We refer the interested reader to p, Q for further details of the data 
and a recent comprehensive theoretical analysis of the same in 0, in which a number of 
non-SM Ansatze are put forward to explain the present experimental anomalies. In the 
meanwhile, the value of the top quark mass has come down marginally; the present world 
average rrit = 175 ± 9.0 GeV has slightly eased the situation for the SM. 

At this particular junction of experiments and SM, it is a fair question to ask: Quo 
Vadis SM? In our view, it is perhaps still too premature to argue persuasively that the 
SM and experiments have parted with each other. This is a tenable point of view as the 
case against the SM is at best circumstantial and by no means impeccable. It should be 
stressed that a good fraction of the LEP data remains to be analyzed. Likewise, lot more 
data will be collected at the SLC collider. The experimental jury is, therefore, still out on 
Rb and Rc and the final verdict on the validity of the SM is yet to be spoken. Improved 
measurements of mw at LEP200, projected to yield a precision of Amw < 50 MeV 0, 
and a projected precision of Am^ = ±4 GeV on the top quark mass at the Tevatron 
collider p , likewise, will also have a direct bearing on the issues being discussed. 

We shall concentrate in these lectures on the physics of the charged current processes in 
B decays, which is not effected directly by the possible modifications of J^*" . However, the 
additional NC couplings, if present, may lead to modifications of some FCNC rare B decay 
rates and distributions, such as in i? — >■ Xi^i~ and B Xuu, in which such couplings do 
play a role. Since, in absolute terms, the possible deviation of the experimental widths in 
question from their SM value is small, we do not expect that such possible modifications 
in J^*-^ will significantly alter the SM-based profile of FCNC B decays. The implications 
of a modified J^*" in B decays must be investigated quantitatively, in case data force 
such a modification. In the CC sector itself, despite suggestions to the contrary 0], we 
shall argue in these lectures that there is no evidence for new physics so far, though this 
may change as data and theory of weak decays get precise and compelling. 

We now proceed to discuss the charged current, J^'" , which is to be derived from the 
£(/, W, B) part of the SM Lagrangian using the mass eigenstates. The CC couplings in 
the SM involve only the left-handed fermions and Z}^. Concentrating on the hadronic 
(quark) part of one can now do the following manipulations: 

V 2 sm Ow i=i 

3 



sin 9w 



^M*^V2PV77'^Vr+V/°"°V/°""^4 + h.c. 



i=l 



tThis is defined as P = (O(expt) - 0(th))/(T(0)(exp). 
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Thus, the charged current J^*", which couples to the W^, is 



^?^= ■ . (^,c,^')l7m^CKM I 5 I , (17) 
singly I ^ 

where we have again dropped the superscript and V^km = V^'^V^"'^^'^ is a (3 x 3) unitary 
matrix in flavour space, first written down by Kobayashi and Maskawa in 1973 [0. It 
is a generalization of the Cabibbo rotation |Tl| for the three-quark-flavour (m, d, s) case. 



invented to keep the universality of weak interactions, which took the form of a (2 x 2) 
matrix by the inclusion of c-quark with the GIM construction |T^, and is called the 



Cabibbo-Kobayashi-Maskawa (CKM) matrix. The charged current Lagrangian has the 
property that it has a. {V — A) structure, hence it violates P and C maximally, conserves 
the electric charge and the lepton- and baryon-number separately, but otherwise there are 
no restrictions on it except that V^j^j^Vckm = 1- In general, violates CP due to the 
possibility of a non-trivial phase in Vckm- 

From this discussion it is clear that the flavour-changing transitions in the SM are 
confined to the CC sector. They emerge in the process of diagonalization of the quark 
mass matrices after spontaneous symmetry breaking, which in turn depend on the Higgs- 
fermion Yukawa couplings. As stated already, the Yukawa couplings in the standard 
model are arbitrary complex numbers. Their adhoc nature in the standard model is in 
all likelihood pointing to the physics beyond the SM, which may help in understanding 
the observed pattern of the fermion masses and mixing angles. How such extensions will 
look like is, however, neither obvious nor the subject of these lectures. Here, we have 
a restricted mandate, namely we will investigate the question if the SM, which comes 
together with these masses and mixings, is a consistent and complete description of data 
or not. 

The matrix elements of V^km determined by the charged current coupling to the 
bosons. Symbolically this matrix can be written as: 



^CKM 



Of the nine elements of VckM; the four involving the u, d, c, s quarks and present in the 
upper left 2x2 block, are studied in decays which are not our principal concern here. 
Their current values can be seen in the PDG review, where also references to the original 
literature can be found 0. 

We shall concentrate here on the remaining five matrix elements in which quarks in the 
third family are involved. Two of the matrix elements involving the h quark, Vub and Vch-, 
have been measured in direct decays of the B hadrons in experiments at the e^e~ storage 
rings CESR, DORIS and LEP. The theory and phenomenology behind their determination 





Vus 


Vub 






Vcb 


Vtd 


Vts 


Vtb 
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will be discussed at some length here. The remaining three elements Vtd, Vts and Vtb can, 
in principle, be directly measured in the decays of the top quark t bW~^ , t sW~^ 
and t — >■ dW~^, respectively. First measurements of Vtb have been reported by the CDF 
collaboration |]13|, through the measurement of the ratio Rtb, 

which is consistent with unity but within experimental errors also consistent with a value 
of Vtb which is considerably less than that. Apart from establishing the dominance of Vtb 
(over the other two matrix elements) by improved measurements of the ratio Rtb | it will 
be difficult in the foreseeable future to get quantitative information on Vtd and Vts from 
decays of top quarks, both due to the scarcity of data involving top quark production 
and, more importantly, the issues having to do with efficient tagging of light-quark jets. 

Fortunately, the matrix elements Vti are also accessible in B and K decays through 
virtual transitions involving the couplings Wtb, Wts and Wtd. Examples of these transi- 
tions in the B system are the \AB\ = 2 processes, B^ - B^ mixings, and |Ai?| = 1 FCNC 
processes, such as rare B decays 6 — > (s, d) + 7 and b — > {s,d) + . We discuss the 
theory and phenomenology of these processes here. Precision experiments involving B de- 
cays and mixings will completely determine the matrix Vq^m, including the CP- violating 
phase, establishing either that the SM provides a consistent theoretical framework for 
describing flavour physics or else that the charged current J"^ must be modified. 

1.1 Some popular representations of the CKM Matrix 

We discuss a couple of popular representations of the CKM matrix. The original parametriza- 
tion due to Kobayashi and Maskawa was constructed from the rotation matrices in 
the flavour space involving the angles 9i {i = 1,2,3) and a phase 6, 

Vkm = R2M. S)Ru{ei, 0)i?23(^2, 0), (20) 

where < < 7r/2, < 6 < 27r, and Rij{6,(j)) denotes a unitary rotation in the 
plane by the angle 6 and the phase 0. The resulting representation is: 



^ Ci -S1C3 -S1S3 

S1C2 C1C2C3 - 82836'^ C1C2S3 + 82036'^ I , (21) 
V S1S2 C1S2C3 + €2836'^ C1S2S3 - €2036'^ 



with Ci = cos 9i, Si = sin 9i. This reduces to the usual Cabibbo form for 62 = O3 = 0, 
with the angle ^1, identified (up to a sign) with the Cabibbo angle. In the PDG review 
0, however, another parametrization is advocated which differs from Vj^M assigning 
the complex phases (dominantly) to the (1,3) and (3,1) matrix elements of V^km- We 
shall not write the PDG representation of the CKM matrix but give instead the simpler 



§A precision of \Vtb\ — 1.0±12% is, for example, projected at the Fermilab Tevatron with an integrated 
luminosity of 2(/6)^^, expected to be collected at the turn of this century S. 
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approximate form, which follows from the a posteriori realization that C13 ~ 1 to a very 
high accuracy, due to the measurement of \Vub\ = S13 ~ 0.003-0.006 (see below). 



' C12C13 S12C13 si3e "^13 

~'Sl2C23 C12C23 S23C13 I . (22) 

\ S\2S2Z — Ci2C23Si3e*'^^^ " C12S23 C23C13 



The three angles called i ^ j-, can be chosen to lie in the range < < 7r/2, making 
all the sines and cosines Cij and Sij real and positive. The KM phase, called for obvious 
reasons ^13, lies in the range < 5i3 < 27r. In the limit ^13 = 6'23 = 0, the third generation 
decouples; identifying the Cabibbo angle with ^12, one recovers the Cabibbo form. To an 
excellent accuracy, which we will discuss, one could set C13 = 1 and C23 = 1. In that case, 
for the PDG parametrization, 

\Vus\ = \si2Cn\ ~ |si2|; \Vub\ = \si3\; \Vcb\ = IS23C13I ^ |s23|- (23) 

These can then be taken as three independent parameters, measured directly in decays, 
together with the phase 613. 

Another approximate but very useful form of the matrix V^km is due to Wolfenstein 
I^, who made the observation that (empirically) the following pattern for the V^km 



matrix elements is suggested by data: 

\Vii\ ~ 1, i = 1...3, 

\Vi2\ ^ IV21I ~ A, 

1^231 - |V^32| ~ 

iV^isI ~ A^, 1^231- A^, (24) 

with A = sin^^c — 0.221. Thus, it is useful to write a perturbative form (in A) for Vqkm- 
Denoting this by Vwolfenstein' 




^Wolfenstein = "A 1 - AX \ . (25) 



Like the previous representations, V^Yolfenstein has also three real parameters called A, X 
and p, and a phase rj. Since we shall be making extensive use of this parametrization, we 
write some relations involving the matrix elements of interest in this representation: 

= X^p^ + v', |^ = Av/(l-p)2 + r/2, (26) 

I ^cb I I ^cb I 

^ IT/ . I 

= 1 , (27) 



\Vtd\ _ j ji-py + v^ \Vt 

\Vub\ V p2 + r/2 ' IK 



cb\ 



and the dominant phases are: 



^{Vub) = Q{Vtd) = -AX't]. (28) 
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It should be recalled that the Wolfenstein parameterization given in Eq. ( P^D is an approx- 
imation and in certain situations in the future it may become mandatory to specify the 
matrix by taking into account the dropped terms in O(A^) in V^volfenstein- ^^e present 
experimental and theoretical accuracy, the representation (^) is entirely adequate and 
we shall restrict ourselves to this form. Further discussions on this point and suggestions 
on improved treatment to include higher order terms in A can be seen in ||15||. 



The four CKM parameters and the six quark masses together with the masses of the 
three charged leptons make thirteen of the nineteen parameters of the standard model; 
the remaining six can be taken as the three gauge coupling constants, the Higgs boson 
mass, the mass of the W boson, and the parameters O^aciQCD) and O^aciEW), which are 
related to the instanton sectors of QCD and the electroweak theory, respectively. Being 
fundamental constants of nature, it is utmost important to measure them as precisely as 
possible and the main goal of flavour physics is to pin down at least the first thirteen of 
them. 



1.2 The CKM unitarity triangles 

The CKM matrix elements obey unitarity constraints, which state that any pair of rows, 
or any pair of columns, of the CKM matrix are orthogonal. This leads to six orthogonality 
conditions. The three involving the orthogonality of columns are listed below, with the 
quark pair in the parenthesis {jk) representing the product of the j'th and fc'th columns: 

ids): Y.y^dV:, = 0, 

i 

(sb): Y.y^sV*, = 0, (29) 

i 

(db): Y.V^<iVi = 0- 

i 

Similarly, there are three more such orthogonality conditions on the rows: 

(uc): = 0, 

j 

(ct): Y^Vc^V* = 0, (30) 

j 

(ut): Y.^u,V* = 0. 
j 

The six orthogonality conditions can be depicted as six triangles in the complex plane 



of the CKM parameter space |T6[. The FCNC transitions in which the pair of quarks 
depicted participate test these triangular constraints directly. Thus, the triangle labeled 
(ds) represents the unitarity constraints on the transition s ^ d, which, for example, one 
encounters in the K^-K'^ transition and in rare K decays such as n^i^i' . The two 

others in this group, namely [sb) and [db), are encountered in the FCNC transitions in 
B decays, such as particle-antiparticle mixings - B^ and B^ - B^ and rare B decays 
B — >■ K* + 7 and B ^ u + respectively. The unitarity conditions on the rows will be 
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Figure 1: The unitarity triangle. The angles a, (3 and 7 can be measured via CP violation 
in the B system and the sides from the CC- and FCNC-induced B decays. 

difficult to test in FCNC transitions shown, as the rates of such transitions are enormously 
suppressed in the SM, and we shall not discuss this set in these lectures. 

The constraint stemming from the orthogonality condition on the ffist and third row 

Vudv:, + v^,vi + v^,v:, = (31) 

has received considerable attention. Since Vud — Vtb — ^ and V^g ~ — Kb, the unitarity 
relation (^) simplifies: 

Vub + = VusVch. (32) 

which can be conveniently depicted as a triangle relation in the complex plane, as shown 
in Fig. |I[ Thus, knowing the sides of the CKM-unitarity triangle, the three angles of the 
triangle a, /? and 7 are determined. These angles are all related to the Kobayashi-Maskawa 
phase 5 (equivalently the phase 5i3 in Vpdg or the phase 77 in Vwoifenstein); they can, 
in principle, be independently measured in various CP-violating B decays. As we shall 
discuss below, the matrix elements Vch and Vub are already known from the CC B decays. 
With more data from B decays and an improved theory one would be able to determine 
them rather precisely. The matrix element Vtd can, in principle, be determined from the 
rare decays b ^ d + •y, b ^ d+ b ^ d + w (and some selected exclusive decays), 

and B^-B^ mixing, which already provides a ffist measurement of \Vtd\ which is, however, 
not very precise. This set of experiments then provides another way of determining the 
triangle, namely by measuring its sides. The unitarity triangle represents an important 
testing ground for the SM flavour physics, in particular in B and K decays. 

It has been pointed out by Jarlskog |]T^, that there exists a large number of different 
parametrizations of the CKM matrix. However, since the phases of the quark fields are 
unphysical quantities, the different parametrizations, emerging from specific choices of 
these phases, must all be equivalent. The parametrization independent quantities are 
the absolute values of the matrix elements \Vij\ (hence also the angles of the unitarity 
triangles) and the area of the unitarity triangles, which is the same for all six triangles 
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and is an invariant measure of CP violation. This can be expressed as 



area[A(CKM)] = ^Sis^ss^ia sin [PDG], 

= IrAX^rj [Wolfenstein], (33) 

for the PDG and Wolfenstein parametrizations of V^km- The Jarlskog invariant denoted 
by the symbol J [0 is twice this area, which in the standard model is typically of O(10~^). 



It is being debated if the intrinsic smallness of J in the standard model is a serious problem 
in explaining the measured baryon asymmetry of the universe (BAU), whose quantitative 
measure is the ratio of the baryon number density to entropy density, 

A. = (34) 

s 

and its present value is = (4 — 6) x 10^^^ 0. Electroweak baryogenesis is a subject 
of great theoretical interest and one which is not quite under quantitative control as 
it requires a deeper understanding of the dynamics of the electroweak theory at high 
energies and temperatures. We refer to a recent review on this subject by Rubakov 
and Shaposhnikov [IT3, where several of the key theoretical concepts and calculational 



strategies are discussed. There is one input to the electroweak baryogenesis issue that can 
possibly come from B decays, namely these decays will probe if nature admits more than 
one mechanism of CP violation. Along the same lines, new sources of CP violation, i.e., 
additional phases, such as are present in multi-Higgs models, may uncover themselves by 
inducing electric dipole moments for the neutron and the electron at a level just above 
their present upper bounds 0, which for sure will not be accommodated by the phase in 
the CKM matrix. 



2 Dominant B Decays in the Standard Model 

With this introduction, we now turn to study questions which are specific to B physics. 
We start by discussing the dominant decay rates which determine the lifetimes of the B 
hadrons, r^, their semileptonic branching ratios Bsl and the charm quark multiplicity in 
B decays (n^), a quantity which has become an important ingredient in understanding 
the semileptonic branching ratio in the standard model. 

The effective lowest-order weak interaction Hamiltonian can be expressed in terms of 
J?,'^ , introduced earlier. 



nw = % (j^^^J^^^^ + h.c.) , (35) 



where Gi? is the Fermi coupling constant. The calculational framework that will be used 
is QCD and we concentrate first on perturbative QCD improvements of the decay rates 
and distributions in B decays. The leading order (in a^) perturbative QCD improvements 
using Tfyy have been worked out in semileptonic processes in [0] - |^5|], which are mod- 



eled on the electromagnetic radiative corrections in the decay of the /x-lepton |]2B[. This 
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argument rests on the property of asymptotic freedom enabling one to use the parameter 
«s(^fe)/7r to do a perturbation expansion of the decay amphtudes. There still remains the 
question of how to relate the perturbative calculations to the actual decays of B hadrons. 
This involves some kind of parton-hadron duality that will be discussed later. For the 
non-leptonic decays, such perturbative corrections are calculated using the renormaliza- 
tion group techniques |]2^ — [Q. We shall discuss some of these techniques here. The 
underlying theoretical framework and its numerous applications in weak decays of the K 
and B mesons have been recently reviewed in a comprehensive paper by Buchalla, Buras 
and Lautenbacher to which we refer for details. 

The quantitative description of the physical decay processes, with hadrons in the ini- 
tial and (in most cases of interest) also final states, requires the knowledge of the wave 
functions (using the quark-parton model language) or hadronic form factors and structure 
functions. The primary task of the theory therefore is to evaluate these non-perturbative 
functions which depend on the strong interaction dynamics. The methods that we shall 
be using here in studying weak decays are based on the lattice QCD framework, QCD 
sum rules, and the heavy quark limit of QCD which allows one to do a systematic ex- 
pansion of decay amplitudes in l/mq, where mg ^ Aqcd, and Aqcd is the QCD scale 
parameter which is typically of 0(200 MeV) 0. The lattice QCD framework aims at 
calculating Green's functions and their S-matrix elements from first principles, i.e. QCD. 
However, in practice, predictions are hampered by the limitations having to do with the 
computing power and/or appropriate lattice formulations. Unquenching lattice QCD is 
one of the foremost technical problems. In addition, simulating B systems directly on the 
lattice with present day technology introduces finite-size effects. In view of this, present 
lattice-QCD estimates in B systems are to be taken with some caution. For results on 
topics relevant for these lectures, we refer to recent reviews ^ |3^ from where we 



shall be drawing heavily. The approach, involving the QCD sum rules [35 , allows one 



to make predictions about current correlation functions which are calculated using the 
operator product expansion (OPE). The results allow themselves to be expressed in terms 
of a limited number of non-perturbative parameters. To extract physical quantities, the 
notion of quark-hadron duality is invoked which enables one to compare suitably weighted 
quantities. While the method is theoretically well-founded, reliable calculations can only 
be made if the higher twist- and higher order QCD corrections have been calculated. This 
is not a limitation in principle but again in practice a trustworthy theoretical treatment 
is available only in a limited number of cases. For a recent review of the applications of 
the QCD sum rules in B decays and discussions of some of the inherent uncertainties, see. 



for example, the review by Braun [36 



Recently, remarkable progress has been made in the formulation of QCD as an effective 
theory in the heavy quark limit, in which the resulting theory shows symmetry properties 
not present in the original QCD Lagrangian. These symmetries enable one to make model- 
independent predictions for hadronic transition form factors involving some exclusive 
B (D, D*, ...) decays at kinematic points where such symmetry relations hold. Some of 



the pioneering work in this direction can be seen in p7|, p8|. In particular, these methods 



have enabled us to determine \Vcb\ with controlled theoretical errors and we shall discuss 
this application here. Along the same lines, equally interesting are the related techniques 
which involve a systematic expansion of the inclusive decay amplitudes in the inverse 
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heavy quark mass |^ - |Q . A very satisfying feature of this framework is that the parton 
model for heavy quark decays emerges as the leading term in a systematic expansion of the 
decay amplitudes. Power corrections to some of the inclusive and exclusive decay processes 
have been calculated in terms of a limited number of non-perturbative parameters. These 
methods, like the QCD sum rules, use operator product expansions which are well defined 
in Euclidean space. To make predictions in time-like regions, some notion of quark-hadron 
duality is again required. Theoretically, this aspect of the effective theories is not quite 
understood. We shall not enter into the technical details surrounding these developments 
and refer to the original literature and some excellent reviews on this subject |4| - g^. 



However, we shall discuss several illustrative applications of this method in these lectures. 

2.1 Inclusive semileptonic decay rates of the B hadrons 

With this theoretical prelude, we start with the assumption that the inclusive decays of 
B hadrons can be modeled on the QCD-improved quark model decays. More specifically, 
while calculating rates, we shall be equating the partial and total decay rates of the B 
hadrons to the corresponding expressions obtained in the parton model, relying on the 



heavy quark expansion ||39[ - 



T{B ^X) =T{b^x) + 0{l/mi) , (36) 

For b quark semileptonic decays involving CC interactions, one has two partonic transi- 
tions: 

b — > ci-ue, (37) 
— > ui^ue. 

There exists a close analogy between the b quark decays and /i decay, fi^ — > e^Uei^^i, 
with the identification: 

[b,{c,u),uej'] ^ [n',e',i?e,i^^]- (38) 

This analogy holds in general; in particular, at the one loop level 0{a) QED corrections 
to /i~ decay and 0{as) QCD corrections to b semileptonic decays are related by simply 
replacing pi, |20|, m 



a 



1 ^ 4 

■ -asTrY,K\ = -Ois, (39) 
"J i=i "J 

where Aj are the Gell-Mann SU{3) matrices, and as is the lowest order QCD effective 
coupling constant, 

«s = — -2 , (40) 

(33-2n^)ln(^^^^) 
^ QCD 

where Uf is the number of effective quarks. The semileptonic decay rates can then be 
read off the expression for the 0{a) radiatively corrected /i-decay rate |^|. The rates for 
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b — > {u, c)ih'e decays, setting me = = 0, are given by the expression: 



■ 2 as{ml) 
1 - o 9{ri 

6 TT 



with To being the normahzation factor in the lowest-order rate 



To 



192^ 



m 



by 



Ti = rrii/mb {i = u, b), and 



fir) 



1 - 8r^ + 8r^ - - 24r^ In r. 



(41) 



(42) 



(43) 



The function g(r) has the normahzation g{0) = vr — ^, and numericaUy (7(0.3) ~ 2.51 



relevant for the b u and b ^ c transitions, respectively |19| 



With A 



QCD 



200 



MeV and Uf = 5, this gives about (15)% corrections to the semileptonic decay widths 
involving £ = e,/i, reducing Tsl compared to the lowest order result vfl = ro/(r). The 
corresponding decrease in the decay width for the semileptonic decay b crur is obtained 
by an expression very similar to the above one in which the r-mass effects are included 
in the phase space and in the QCD corrections. 



r(6 CTUr) = ToP{Xc, Xr, 0) 



1 + '^^^^^\ {Xc,Xr,0) 
07T 



(44) 



where P{xi,X2,X3) is the well known three-body phase space factor given for arbitrary 
masses Xj = mi/rrih by 



P{Xi,X2,X3) =12 J 



(X2+X3)2 



ds ^ 
s 



s-xl- xl){l + xl- s)w{s, xl, xl)w{s, xl, 1) , (45) 



w{a, b, c) = (a^ + b^ + - 2ab - 2ac - 26c) 



1/2 



(46) 



The function g{xi,X2, x^) has been calculated for arbitrary arguments in p3[ in terms of a 
one-dimensional integral. The functions P{xi, 0, 0) and g{xc, 0, 0) go over to the functions 
/(r) and {—)g{r), respectively, given above for the massless lepton case. The numerical 
values for g{xc, Xr, 0) and g{xc, 0, 0) are tabulated in [^. For the default value Xc = 0.3, 
one has g{xc, Xr, 0) = —2.08, yielding about a 12 % decrease in T{b — > cri/^) compared to 



Tsl{b CTUr) as a result of the leading order QCD corrections p3[. For more modern 
calculations of the decay rate Tsiib — >• crur), see ||2 



2.2 Inclusive non-leptonic decay rates of the B hadrons 

The dominant CC-induced non-leptonic and semileptonic decays of B hadrons are gov- 
erned by the effective Lagrangian, 

C^ff = -4^K:,K, [Ci(/i)Oi(/i) + C2i^i)02{fi)] 
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+ h.c. 



(47) 



and we have just discussed the 0{as) renormahzation effects to the matrix elements of 
the semileptonic piece in C Here Oi and O2 are the colour-octet and colour-singlet 
four- Fermi operators, respectively (here a and (3 are colour indices), 



O2 = {daUa)L{cpbp)L, 



(45 



and qi = 1/2(1 — 75) denotes a left-handed quark field. The operators are related to 
the corresponding fields Oi by the relacement d ^ s. The octet-octet (Ci) and singlet- 
singlet {O2) operators emerge due to a single gluon exchange between the weak current 
lines (quark fields) and follow from the colour charge matrix (T^j) algebra: 



^ik-'-jl 



(49) 



Here, A^c = 3 for QCD. The Wilson coefficients Cj(/i) are calculated at the scale = mw 
and then scaled down to the scale typical for B decays, n = 0{mf,), using the renormal- 
ization group equations, which brings to the fore the influence of strong interactions on 
the dynamics of weak non-leptonic decays. Without QCD corrections, the two Wilson co- 
efficients have the values Ci(mvy) = 0, C2(mvi/) = 1- Since the operators Oi and O2 mix 
under QCD renormahzation, it is convenient to introduce the operators 0± = {02^0i)/2 
having the Wilson coefficients C± which renormalize multiplicatively |2^. The results are 



now known to two-loop accuracy pO[]: 



(0) 



1 + 



4vr 2/?o 



7^^ 
7?^ 



as{mw) ^ 



(50) 



where the multiplicative factor in this expression represents the solution of the RG equa- 
tions in the leading order QCD p7| . 



as(M, 



w. 



d± 



(51) 



and the exponents have the values d+ = 7+^(2/^0)5 d- = '-/^^^ /{2/3o). The quantities 7+"* 
are the coefficients of the anomalous dimensions involving the operators 0± (and 0±), 



with 



(0) A (0) 

7; -4, f_' 



(1) '7,4 (1) 

=-7+g^/, 7- 



(52) 
(53) 
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in the naive dimensional regularization (NDR) scheme, i.e., with anticommuting 75. The 
j3i are the first two coefficients of the QCD /5-function, and they have the values 



f3o = ll-\nf, /5i = 102-yny. (54) 

Finally, the functions B± are the matching conditions obtained by demanding the equality 
of the matrix elements of the effective Lagrangian calculated at the scale /i = my/ and in 
the full theory (i.e., SM) up to terms of 0(as(m^)). They have the values: 

B. ^ ±B^, (55) 
The constant B and the two-loop anomalous dimension 7^'' are both regularization- 



scheme dependent. In the NDR scheme one has B = 11. Following [^, we define a 
scheme-independent quantity -R±, 

^(0) .-,(1) o . 



V^f A 

in terms of which the Wilson coefficients read 



C±(/i) = L±(/i) 



^ + + -^^^ 



(57) 



In this form all the scheme-dependence resides in the coefficients B± which is to be can- 
celled by the scheme-dependence of the matrix elements of the corresponding operators. 

In addition to the decays 6 ^ c + ud, b c + us and 6 — > c + iu^, which are described 
by the effective Lagrangian (^Tf ), there are other decays involving the CC transitions 
b —>■ uX, b —>■ (c, u) + CS and b —>■ (c, u) +cd, which are not included in this Lagrangian. In 
a systematic treatment involving QCD renormalization, one has to enlarge the operator 
basis to include these transitions and the so-called penguin operators. We shall return to 
a discussion of this part of the Lagrangian later in these lectures as we discuss rare B- 
decays, where the operator basis will be enlarged and the corresponding Wilson coefficients 
calculated in the leading logarithmic approximation. 

We now discuss the semileptonic branching ratio SgL for the B mesons and to be 
specific will consider the case i = e, fi. This branching ratio is to a large extent free of 
the CKM matrix element uncertainties but requires a QCD-improved calculation of the 
inclusive decay rates, Fg^ , discussed above, and F^q^^, 



with 



tot 



(B) = T{B^Xiiye)+T{B^X,X) + T{B^X,,X) 



-e,fi,T 



+ T{B XuX) + F(5) (Penguins) . (59) 
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In the spirit of the parton model, we shall equate r{B XcX) = r{b — > cud) + r{b 
cus), noting that the so-called VT-annihilation and VT-exchange two-body decays are 
expected to be small in inclusive B decays. This will be quantified later as we discuss the 
lifetime differences among B hadrons which arise from the matrix elements of the operators 
representing these contributions. The corrections for the decay widths r(6 — cud) and 
r(6 — cus) are identical neglecting rriu and m^, and so their contributions can be described 
by similar functions. The resulting next-to-leading order QCD corrected sum can be 
expressed as: 



r(6 ^ cud) + T{b 



cus) = TqP{xc, 0, 0) 

X 



27r 



+ 



(2L(/i)^i?+ + L(/.)ii?_) 
3 



1 



+ ^{L{fi)l - L(^)^_)(ci2(x„/^) - 121n-^) 



rrib 



3Tor]{fi)J{xc,fi) 



with ri{fi) representing the leading order QCD corrections. The scheme independent R± 
come from the NLO renormalization group evolution and are given by |^ 



R. 



10863 - 1278n/ + 80nj 
6(33 - 2nj)2 ' 
15021 - 1530n/ + 80n} 
3(33-2n/)2 



(61) 



For Uf = 5, i?+ = 6473/3174, i?_ = -9371/1587. Note that the leading dependence of 
L{fi)± on the renormalization scale fj, is canceled to 0{as) by the explicit /^-dependence 
in the a^-correction terms. Virtual gluon and Bremsstrahlung corrections to the matrix 
elements of four fermion operators are contained in the mass dependent functions Cij{x). 
The analytic expressions for the functions Cii(x), Ci2{x), 022(2;) are given in where also 
their numerical values are tabulated. For our default value Xc = 0.3, these coefficients 
have the values: 



cii(0.3)~2 ci2(0.3,/i = mb) ~ -10 C22(0.3) ~ -1 . 



(62) 



Lumping together all the perturbative and finite charm quark corrections in a multiplica- 
tive factor Ac{mb, Xc, as{mz)), the perturbatively corrected decay width can be expressed 
as: 

r(6 — > cud) + r(6 cus) = 3roP(a;c, 0, 0) [1 + Ac(xc, mi,, as{mz))] ■ (63) 

For the central values of the parameters used here {mf, = 4.8 GeV, Xc = 0.3, fj, = nib and 

as(mz) 



0.117), the QCD corrections lead to an enhancement ^1 
Ac(mfe, Xc, as{mz)) = 0.17. 



(64) 



(60) 
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Out of this, the bulk is contributed by the leading log factor 



r/(/i) - 1 = - (24 + Ll 



1 = 0.10 



(65) 



Next, we equate V{B Xcc) = r(6 — ccs) + r(6 — > ccd) and discuss the perturbative 
QCD corrections to the decay width V{b — >• ccs) and V{h — >• ccd). Neglecting rrid and m^, 
an assumption which has been found to be valid to a high accuracy in WM, the corrections 
in the two decay widths are identical and the result can be written in close analogy with 
the ones for the decay widths V{h ^ cus) discussed above. 



V{h ^ ccs) + r(6 — i> ccd) 



X 



+ 



2L(^)^ + L(^)^ 



as{Mw) - as (/i) 
27r 



(2L(^)^i?+ + L(/i)2 i?_ 



37r 



3 3 

-(L(/i)+ - L{ii)J)'^kii{xc, /i) + -{L+ + LJ)'^k22{xc 



(66) 



The functions kij{mf,, Xc, ctsimz)) have been calculated and their numerical values are 
tabulated in |Q. For = 4.8 GeV, Xc = 0.3 and as{mz) = 0.117 and fi = rrih, they 
have the values: 



kii{0.3) = 6.44 A;i2(0.3, fi = rrib) = 0.82 A:22(0.3) = 2.99 



(67) 



showing that these coefficients are rather sensitive functions of Xc- The corresponding 
numbers for kij for the choice ms/rrih = 0.04 are given in |Q. Again, lumping to- 



gether all the perturbative and finite charm quark corrections in a multiplicative factor 
Acc{rni„ Xc, asimz)) , the perturbatively corrected decay width can be expressed as: 



T{b ccs) = 3ToP{xc, Xc, Xg) [1 + Acc(xc, m^, as{mz))] . 
For the values of the parameters used here (m^ = 4.8 GeV, Xc = 0.3 and as{mz) 



the QCD corrections lead to an enhancement p3 

Acc{mb,Xc,as{mz)) = 



0.37. 



(68) 
0.117), 

(69) 



Out of this, the bulk is contributed by the next to leading order correction. This is by far 
the largest correction to the inclusive rates we have discussed so far. Using pole quark 
masses and the renormalization scale 



rrib, one gets 
ccs){NLO) 



r{b^ccs){LO) 



1.32 ±0.07 



(70) 



The NLO corrections go in the right direction in bringing theoretical estimates closer to 
the experimental value for the semileptonic branching ratio. However, this will also lead to 
enhanced charmed quark multiplicity (uc) in B decays. We shall return to a quantitative 
analysis of SgL (-B) and (hc) at the end of this section. 

We now turn to a discussion of the CKM-suppressed and penguin transitions con- 
tributing at a smaller rate to r^Q^{B). They are of two kinds: 
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• r{B Xu + X), which is suppressed due to the CKM matrix element \Vub\, with 
the rate depending on |V^6p, and 

• r (5) (Penguin): The so-called penguin transitions b ^ s + X , where X = cc and 
^ = 9 (QCD penguins), X = 7 (electromagnetic penguins), X = i^(i~,ui) (elec- 
troweak penguins). 

There are even smaller transitions involving h ^ d + X, as well as a host of other rare 
decays but we shall neglect them in estimating the total decay width for obvious reasons. 
We list below the numerical contributions where in all entries involving 5 — > m + X 
transitions, we have set = 0.08, corresponding to the present central value of 

this ratio The contribution of 6 — > uud is calculated without penguins and the 

contribution to the h —>■ ccs rate given below is due to interference of the leading current- 
current type transitions with the penguin operators. 

T{b^uJ2 ~ O.OlSro r(6 ucs) ^ O.OlOFo (71) 

I 

r{b uud) ^ 0.022ro ATpenguinib ccs) ^ -0.04iro (72) 

T{b sg) ^ 5.0 X IO^^Fq T{b 57) ^ 8.0 x IO^^Fq (73) 

T{b svv) ^ 1.2 X lO'^Fo T{b stT) ^ 5.0 x IO^^Fq (74) 

Note that the contribution due to the interference with the penguin transitions in 6 — »• ccs 
is negative and it tends to cancel the other small contributions listed above in the total 
non-leptonic width. The sum of all these contributions add up to 

T{B ^ Xu + X) + r (5) (Penguins) - 1.25 x IO-^Fq , (75) 

and hence not of much consequence for the semileptonic branching ratio or the B hadron 
lifetime estimates. Rare B decays are either constrained by direct experimental searches or 
indirectly through the measurement of the branching ratio B{B —>■ Xg + 'j) which excludes 
most of the allowed parameter space where these decay modes may have appreciably 
larger branching ratios. It may be parenthetically added here that the chromomagnetic 
Q CD-penguin contribution b ^ s + g could be significantly large in some extensions 
of the SM, in particular the minimal supersjTiimetric standard model (MSSM) []55| , 56| . 



This is sometime suggested as a possible source of enhanced non-leptonic decay width. 
There exist already an experimental bound on this decay branching ratio 0], obtained 



through the decay chain b—i'S + g^Xg + cpoia the assumption that the 0- energy 
spectrum in inclusive B decays is rather hard, similar to the photon energy spectrum 
in electromagnetic penguin decays B + 7. This quasi two-body i^^-spectrum 

is suggested in ^7j. However, in our opinion, such an assumed i^^-spectrum is rather 
unrealistic as it is unlikely that the fragmentation of an 0(5 GeV) s + g system will give 
rise to dominantly two-body or quasi-two body final states, hence the experimental bound 
so-obtained is not very compelling. 
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2.3 Power corrections in TsiiB) and TntIB 



Before we discuss the numerical results for i3gL, we include the 0{l/m\) power corrections 
in the inclusive partonic decay widths, which have been calculated using the operator 
product expansion techniques [^- [|^] and they constitute the first non-trivial corrections 
to the parton model results. Taking the typical hadronic scale to be 0(1) GeV, one expects 
as a ball-park estimate 0(5)% corrections in the inclusive rates. This is now discussed 
more quantitatively. 

In HQET, the 6-quark field is represented by a four-velo city- dependent field, denoted 
by hv{x). To first order in l/mf,, the 6-quark field in QCD b{x) and the HQET- field by{x) 
are related through: 



b{x) 



l + i 



2mb 



The QCD Lagrangian for the b quark in HQET in this order is: 



^HQET 



b^iv. pby + b^— b^ 



Zhb 



h'Jv- 



af3 



a 



Arrih 



-b, + 



mt 



(76) 



(77) 



where Zb is a renormalization factor, with Zb{fi = rub) = 1 and p = D^'j^, with 
being the covariant derivative. The operator b^{i p)%^/2mb is not renormalized due to 
the symmetries of HQET. (In technical jargon, this is termed as a consequence of the 
reparametrization invariance of C^^^^.) With this Lagrangian, it has been shown in 
PP| - that in the heavy quark expansion in order (1/m^), the hadronic corrections 
can be expressed in terms of two matrix elements 



{*) 



) = 2dB(,)mB(*)\2 



(78) 
(79) 



where F^'^ is the gluonic field strength tensor, and the constants dQ(*) have the value 3 and 
— 1 for B and 5*, respectively. The constant A2 can be related to the hyperfine splitting 
in the B mesons, which gives: 



A, 



B* 



m 



I) 



0.12 GeV^ 



^0) 



The other quantity Ai is just the average kinetic energy of the b quark inside a B meson 
and has been estimated in the QCD sum rule approach yielding Ai = —(0.5 ± 0.1) 
GeV^ [^. As we shall discuss later, this parameter influences the lepton- and photon- 
energy spectrum in the decays B —>■ Xiue and B ^ Xg + 'y, and has also been estimated 
using these spectra |^ Taking into account these corrections, the semileptonic and 
non-leptonic decay rates of a 5 meson B —* Xiui and B X^X can be written as 



T{B^Xdy() = rW/(r,) 



2 a, 

3 TT 




6(1 - r,Y A2 



2m? 



2ml 



fire 



mt 



a. 



m( m'l' 



(81) 
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and 



3r(°) [r/(/i)J(/i) (^1 
-(L+(/i)2-L_(/if 



Ai 
2ml 



+ 



3X2 
2ml 



fire 



mt 



a. 



fire 



mt 



mt mv 



(82) 



where the product ri{fi)J{fi) denotes the NLO corrected resuh for the partonic decay 
discussed above in (|60|) , to which Eq. (^) reduces in the hmit Ai = A2 = 0. 

The decay rates depend on the quark masses, which unhke lepton masses, do not 
appear as poles in the S-matrix nor do the quarks exist as asymptotic states. They 
are parameters of an interacting theory and hence subject to renormahzation effects. 
Consequently, they require a regularization scheme, such as the MS scheme, and a scale, 
where they are normalized, to become well-defined quantities. For example, the quark 
masses in the so-called MS scheme and the pole masses (OS scheme) are related in the 
leading order 



mQ[mQ) = mg 



1-4- 



(Svr) 



+ 



83) 



In HQET, quark masses can be expressed in terms of the heavy meson masses mM 
and the parameters Ai, A2 and a quantity called A, where 



J Xi + dMX2 , 
rriM = mQ + A h ... 



2m, 



Q 



With this the decay rates then depend on the QCD-related parameters, Aqcd, Ai, A2 and 
A. Since quark masses are scale- and (regularization) scheme-dependent quantities, this 
also holds for the parameter A. As a consequence of this, the decay rates become scale- 
and scheme-dependent. This will reflect itself through increased theoretical dispersion 
on various physical quantities. More precise predictions for the decay rates and related 
quantities (such as Bsl and {ric)) require knowledge of higher order corrections, which are 
unfortunately not known. In the context of the heavy quark expansion, one can relate 
the parameters Ai and A2 to the quark and hadron masses: 



nib - mc 



A1 + 3A2, 1 1 , 1 , 

ms-niD^ ( ) + 0[—) 

I mh mr m^ 



^5) 



and the quark mass differences can then be calculated knowing Ai and A2, giving [mt — 
m,) = (3.4 ± 0.03 ± 0.03) GeV This difference, which determines the inclusive 

rates and shape of the lepton energy spectrum in semileptonic decays, has also been 
determined from an analysis of the experimental lepton energy spectrum in B decays, 
yielding [mb — m,) = 3.39 ± 0.01 GeV for the pole masses [0], in excellent agreement 
with the QCD sum rule based estimates. The quark masses themselves have considerable 
uncertainties. We shall not discuss them here. For a comprehensive discussion of the 
quark masses and their current values, we refer to the reviews by Leutwyler |^ and by 
Manohar in the PDG review of particle properties 
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2.4 Numerical estimates of BsiiB) and 



{rir 



The theoretical framework described in the previous section can now be used to predict 
two important quantities in B decays Bsl{B) and {ric), which have been measured. Con- 
cerning Bsl{B), there is some discrepancy between the two set of experiments performed 
at the T(4S') and at the Z° resonance, ahhough it must be stressed that these experi- 
ments measure a different mixture of B hadrons. The present measurements are reviewed 



recently in 54 



Bsl{B) = (10.56 ±0.17 ±0.33)% at T(45) 
BsLiB) = (10.89 ±0.18 ±0.24)% at Z° 
(ric) < 1.16 ±0.05 at T(45). 



^6) 



The first error quoted for Bsl is experimental and the second is an estimate of the model 
dependence. The number for (ric) is obtained by summing over various bound and un- 
bound charmed hadron states in B decays, and the inequahty is to due to the present 
upper limit on the inclusive branching ratio B{B Vc^) < 0.018. 

The theoretical predictions for these quantities at present are somewhat fuzzy due 
to the uncertainties in the input parameters. We shall rely here on a recent theoretical 



update by Bagan et al. ||52| [see also the recent analysis by Neubert and Sachrajda p^ ] 
where the following ranges of parameters have been used: 



m6(pole) = 4.8 ± 0.2 GeV; a,(m^) = 0.117 ± 0.007, mb/2 < fi< 2mb 
Al = -(0.5±0.1) GeV^; As = 0.12 GeV^ 



^7) 



Their analysis leads to the following values 1^2 

Bsl 

and 



12.0 ±0.7 ±0.5 ±0.2+°;^) 



Bsl = (11.3 ±0.6 ±0.7 ±0.2: 



^9) 



(90) 



using the pole (also called on shell OS) and MS masses, respectively. The errors are from 
A(mb), Aas{mz), A(Ai), A(r(_B X^c) and from the scale(= /x) variation, respectively. 
An estimate of the present theoretical uncertainty can be obtained by adding these errors 
in quadrature (a reasonable procedure since most errors given above are independent), 
yielding 



Bsl = (12.0 ±1.4)%, 
Bsl = (11.3 ±1.6)%. 

This shows that this quantity has still significant theoretical spread: 

5Bsl{B) - ±0.15Bsl{B) ~ ±0.018 . 



(91) 



(92) 



The corresponding estimates for {ric) and (ric) in the OS and MS schemes, respectively, 
are [|52l : 



(ric) = (1.24 ±0.05 ±0.01) and (n^) = (1.30 ± 0.03 ± 0.04), 



(93) 
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where the first error in both the schemes is due to Am;,, and the second error combines 
the uncertainties in the rest of the parameters. The estimates (0) and (0) should be 
compared with the experimental measurements given in Eq. (|86D . A comparison of these 
theoretical estimates (eqs. (pQ]), (pOD and (0)) and data on (n^) and Bsl is shown in 
Fig. Ij. It should be pointed out that in the analysis of Bagan et al. for Bsl and 
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Figure 2: The charm content ric vs. Bsl B decays. Solid lines: theoretical predictions 
in the OS scheme for 0.23 < mc/mh < 0.33; dashed lines: the same in the MS scheme for 
0.18 < fhc^fhc) / fhiyijhh) < 0.28. Shaded area: theoretical predictions in the OS scheme by 
varying over the parameters as discussed in the text. The data point is from (Figure 
taken from p^). 



(ric), possible contributions from the so-called spectator effects (ly-annihilation and W- 
exchange) have been neglected. Such effects depend on the wave-functions at the origin, 
using the quark model language, and are proportional to which could lead to some 
0(5)% effects, well within the present theoretical noise 6{Bsl) discussed above. All in all, 
it is fair to conclude that within existing uncertainties, the current theoretical estimates 
for Bsl and do not disagree significantly with the experimental values, though there is 
a tendency in these estimates to yield somewhat larger value for n^. More work is needed 
to make these comparisons precise. While still on the same subject, we note that the 
pattern of power corrections predicted by the short-distance expansion of the heavy quark 
effective theory in inclusive non-leptonic decays has been questioned in It is argued 
on phenomenological grounds that data on b- and c-decays are better accommodated if 
some account is taken of the linear (l/mg) corrections, which are absent in the HQET 
approach. Naive identification of the ttiq factor in the inclusive decay widths with mf^, 
with rriH being the heavy hadron mass, brings data and theoretical estimates in better 
agreement! The presence of l/rriQ corrections could be related to the breakdown of local 
duality in non-leptonic decays. The real question is if one could derive this suggested 
pattern of power corrections in QCD. 
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While the suggestion of AltareUi et al. is interesting and it raises new theoretical 
issues as to the reliability of Tj^l calculated in the HQET approach, their arguments 
are purely phenomenological. It is fair to say that a big step in reducing the present 
theoretical uncertainties will be the completion of all the perturbative NLL corrections to 
the dominant decays discussed above. Parts of them, the so-called gluon bubble graphs 
having an arbitrary number of fermion loops, which take into account the effects of the 
running of on Bsl in all orders in perturbation theory, are available in the literature 
371. This is related to the issue of defining a factorization scale, discussed earlier in a 



more general context by Brodsky et al. [^, which at present is a big uncertainty in this 
approach, as can be seen in Fig. |[ 

2.5 B-Hadron Lifetimes in the Standard Model 

A matter closely related to the semileptonic branching ratios is that of the individual 
B hadron lifetimes. The QCD-improved spectator model gives almost equal lifetimes. 
Power corrections will split the 5-baryon lifetime from those of B^, B^ and Bg. However, 
first estimates of these differences are at the few per cent level The experimental 

situation has been summarized as of summer 1995 in [|69|] : 

4^ = 1.02 ±0.04; 4^ = 1-01 ±0.07; 4^ = 0-76 ± 0.05 . (94) 
r{Bd) T[Bd) T{Bd) 

Since then, a new result for the lifetime of the A;, baryon has been reported by the CDF 
collaboration, T{Ab)/T{Bd) = 0.85 ± 0.10 ± 0.05 [0 , which is still less than the other 
two ratios involving the i?-meson lifetimes but reduces the gap in the Af, and i?-meson 
lifetimes. 

This subject has received renewed theoretical attention lately ^ |72|, in which the 
possibly enhanced roles of the four-Fermion operators between baryonic states has been 
studied. We recall that such operators enter at 0(1/ ml) in the heavy quark expansion 
discussed above 0] . In this order, there are four such operators, which using the notation 
of [|^, can be expressed as: 

7v_A = {bLlt,taqL){qLYtabL), 

TLp = {bRtaqL)iqLtabR), (95) 

where ta are generators of colour SU (3). The matrix elements of these operators between 
various i?-meson and A^-baryons are in general different and this contribution will thus 
split the decay widths of the various B hadrons. In general, the operators ( P5| ) introduce 
eight new parameters corresponding to the matrix elements of these operators. In the 



large- A^c limit, however, it has been argued in |Q that the i?-mesonic matrix elements of 
the operators {Bg\Oy_A\Bg) and {Bg\Og_p\Bg) are the dominant ones. While accurate 
numerical estimates require a precise knowledge of these matrix elements, one expects that 
they give rise typically to the spectator-type effects (using the parton model language): 

^ ^ 5% , (96) 

Ttot 
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with Jb of order 200 MeV. In the case of A;, baryons, one can use the heavy quark spin 
symmetry to derive two relations among the operators between the A;, states. The problem 



is then reduced to the estimate of two matrix elements which in ||6J] are taken to be the 
following: 

and 

{Ab\dl_^\At) = -B{A,\Ol_A\Ab) , (98) 

The operator Oy-A is a linear combination of the operators Ty-A and Oy-A introduced 
earlier, Oy^A = 27y_^ + 3Cv-A) following from colour matrix algebra and r{Ab/B, 



is the ratio of the squares of the wave functions which can be expressed in terms of the 
probability of finding a light quark at the location of a 6 quark inside A;, baryon and the 
B meson, i.e. 



(99) 



One expects i? = 1 in the valence-quark approximation. However, the ratio r{Ai,/ Bq) has 
a large uncertainty on it, ranging from r{Ab/Bq) ^ 0.5 in the non-relativistic quark model 
|73| to r{Ab/Bq) = 1.8 ± 0.5 if one uses the ratio of the spin splittings between Ej, and 
baryons and B and B* mesons, as advocated by Rosner and using the preliminary 
data from DELPHI, m(S^) - m(Sfe) = (56 ± 16) MeV WM. 



Using the ball-park estimates that B and r{Ah/ Bq) are both of order unity yields for the 
lifetime ratio r(Afe)/r(i?^) > 0.9 0|, significantly larger than the present world average. 



Reliable estimates of these constants can be got, in principle, using lattice-QCD and QCD 
sum rules. Very recently, QCD sum rules have been used to estimate (Af,|(!)y_^|Aft) and 
5, yielding (AftlO^.^IA;,) = (0.4 - 1.2) x lO'^ GeV^ and B = 1.0 |7|. This corresponds 
to the parameter r{Ai,/Bq) having a value in the range r{Ai,/Bq) ~ 0.1 — 0.3, much too 
small to explain the observed lifetime difference. One must conclude that the lifetime 
ratio T{Ab)/T{Bd) remains a puzzle. New and improved measurements are needed, which 
we trust will be forthcoming from HERA-B and Tevatron experiments in not-too-distant 
a future. 



2.6 Determination of \Vcb\ and \Vub\ 

The CKM matrix element Vet can be obtained from semileptonic decays of B mesons. 
We shall restrict ourselves to the methods based on HQET to calculate the exclusive 
semileptonic decay rates and use the heavy quark expansion to estimate the inclusive rates. 
Concerning exclusive decays, we recall that in the heavy quark limit {rrib oo), it has 
been observed that all hadronic form factors in the semileptonic decays B —>■ {D, D*)lvi 
can be expressed in terms of a single function, the Isgur- Wise function . It has been 
shown that the HQET-based method works best for B — »• D*lv decays, since these are 



unaffected by l/mg corrections |73, |78[. Since the rate is zero at the kinematic point 
= 1, one uses data for > 1 and an extrapolation procedure (discussed below) to 
determine ^(l)|Vcfe| and the slope of the Isgur- Wise function p^. 
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Using HQET, the differential decay rate in B D*ii'£ is 
dT{B D*lv) 



G'f 
48^ 



mD'')'^m%.r]\Vuj'^ - l(cj + 1)^ 
Auj 1 — 2ujr + 



(100) 



X 1 + 



UJ + l 



where r = ttid* /niB, uj = v ■ v' {v and v' are the four- velocities of the B and D* meson, 
respectively), and is the short- distance correction to the axial vector form factor. In 
the leading logarithmic approximation, this was calculated by Shifman and Voloshin some 
time ago - the so-called hybrid anomalous dimension [|79l. In the absence of any power 
corrections, ^{uj = 1) = 1. The size of the 0(l/m^) and 0(l/m^) corrections to the Isgur- 
Wise function, ^(cu), and partial next-to- leading order corrections to rjA have received 
a great deal of theoretical attention, and the state of the art has been summarized by 
Neubert [^] and Shifman ||4^. Following [^], we take: 



Va 



1 + 5(l/m2) = 
0.965 ±0.020 



0.945 ±0.025 



(101) 



This gives the range [^]: 

= = 0.91 ±0.04. (102) 

Recently, the quantity rjA, and its counterpart for the vector current matrix element renor- 
malization, riy, have been calculated in the complete next-to- leading order by Czarnecki 



1 81, getting 



VA 

Vv 



0.960 ±0.007 
1.022 ±0.004 



(103) 



The NLO central value for 77^ is in agreement with the estimate of the same given in 
eq. ( |101| ) but the error on it is now reduced by a factor of 3. So, the error on J-'{1) is now 
completely dominated by the power corrections in ^(1). 

Since the range of accessible energies in the decay B D*iu is rather small 
(1 < uj < 1.5), the extrapolation to the symmetry point can be done using an expansion 
around u = 1, 

Tito) = J^(l) [1 - p2(cj - 1) ± c{uj - If + ...] . (104) 

It is usual to use a linear form for extrapolation with the slope left as a free parameter, 
as the effect of a curvature term is small [^. The present experimental input from the 
exclusive semileptonic channels is based on the data by CLEO, ALEPH, ARGUS, and 

|83| to which we refer for details and references to 

The statistically 



DELPHI, which is summarized in 



the experimental analysis. For the updated ARGUS numbers, see K 
weighted average used in the analysis ||8^ is: 



iV^fcl ■ J^(l) = 0.0353 ±0.0018 



This agrees with the numbers presented by Skwarnicki last summer [P 



\Vcb\-H^) 



0.0351 ±0.0017 
0.87 ±0.10 . 



(105) 



(106) 
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Using ^^-"(1) from Eq. ( |102D , gives the following value: 



iKfel = 0.0388 ± 0.0019 (expt) ± 0.0017 (th). 



(107) 



Combining the errors linearly gives \Vcb\ = 0.0388 ± 0.0036. This is in good agreement 
with the value \Vcb\ = 0.037lo;oo2 obtained from the exclusive decay B — > D*ii'i, using a 
dispersion relation approach . 

The value of | Kfel obtained from the inclusive semileptonic B decays using heavy quark 
expansion is quite compatible with the above determination of the same. The inclusive 
analysis has the advantage of having very small statistical error. However, as discussed 
previously, there is about 2a discrepancy between the semileptonic branching ratios at 
the T(4S') and in Z° decays. Using an averaged value for the semileptonic decay width 
from these two sets of measurements gives: T{b ciui) = Bsl/ {tb) = (67.3 ± 2.7) ns^^"*^^ 
H, where (tb) = {tb- + tb^>)/2 = 1.60 ± 0.03 ps and {tb) = 1.55 ± 0.02 ps have been 
used for the T(4S') and decays, respectively, and the error has been inflated to take 
into account the disagreement mentioned. This leads to a value ^4 



\Vcb\ = 0.0398 ± 0.0008 (expt) ± 0.004 (th) 



(108) 



where the theoretical error estimate (±10%) has been taken from Neubert |8^. While 
still on the same quantity, it should be noted that Vainshtein has quoted a much smaller 
theoretical error, |^ : 



IK 



cb\ 



0.0408 


'B{B^XJu,y 


1/2 


■1.6ps" 




0.105 




- Tb . 



1/2 



;i.0±0.03(th)) . 



(109) 



For further discussion of these matters we refer to 
values for \Vcb\ and the Wolfenstein parameter A: 



5D|. We shall use the following 



cb\ 



0.0388 ± 0.0036 



A = 0.80 ±0.075 



(110) 



Up to recently, iKfe/Kfel was obtained by looking at the endpoint of the inclusive 
lepton spectrum in semileptonic B decays. Unfortunately, there still exists quite a bit 
of model dependence in the interpretation of the inclusive data by themselves yielding 
iKiftl/lKbl = 0.08 ± 0.03 [^, |4|. A recent new input to this quantity is provided by the 
measurements of the exclusive semileptonic decays B — >• (7r,p)£z/£ [^, The extracted 
branching ratios depend on the model used to correct for the experimental acceptance 
but they do provide some discrimination among the various models. In particular, models 
such as that of Isgur et al. ||88|, which give values well in excess of 3 for the ratio of the 



decay widths T{B^ p-i+u)/T{B^ 



TT 



^^^z/), are disfavoured by the CLEO data which 



yield typically 1.7 ± 1.0 for the same quantity, with some marginal model dependence. 



The models with larger values for the ratio T{B 



-u)/T{B° TT-i+u) also lead to 



softer lepton energy spectrum in inclusive B decays, requiring larger values of iKbl/lKfel. 
A recent estimate of this ratio in the light-cone QCD sum rule approach gives 2.4 ± 
0.8 and yields a branching ratio B{B — > piuf) = (19 ± 7)|Vufep, giving a value of 
(1.95±0.72) X 10~^, for \Vub\ = 0.08, in agreement with the experimental numbers, which, 
however, have significant errors. Excluding models which fail to reproduce the exclusive 
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data from further consideration, measurements in both the inclusive and exclusive modes 
are compatible with 

=0.08 ±0.02, (111) 



cb 



which is also the value adopted by the PDG Review 0. This gives a constraint on the 
Wolfenstein parameters p and 77: 



p2 + ^2 = 0.36±0.08 . (112) 

With the measurements of the form factors in semileptonic decays B (vr, p, uj)lv^, one 
should be able to further constrain the models, thereby reducing the present theoretical 
uncertainty on this quantity. 

We summarize this section by observing that the bulk properties of B decays are largely 
accounted for in the standard model. On the theoretical front, parton model estimates 
of the earlier epoch have been replaced by theoretically better founded calculations with 
controlled errors, though this point of view has not found universal acceptance [^! In 
particular, methods based on HQET and heavy quark expansion have led to a quantitative 
determination of \Vch\ at ±10% accuracy, which makes it after \Vud\ and |K(s|, the third 
best measured CKM matrix element. The matrix element \Vuh\ has still large uncertainties 
(±25%) and there is every need to reduce this, as it is one of the principal handicaps at 
present in testing the unitarity of the CKM matrix precisely (more on this later). The 
quantities Bsl-, (nc), and the individual 5-hadron lifetimes are not in perfect agreement 
with data but the present theoretical estimates have still large uncertainties to abandon 
the SM. A completely quantitative comparison requires the missing NLL corrections and 
in the case of lifetime differences better evaluations of the matrix elements of four-quark 
operators, which we hope will be forthcoming. Finally, we stress that it will be very 
helpful to measure the semileptonic branching ratios Bsl for the baryons. With the 
lifetimes of the B hadrons now well measured, such a measurement would allow to compare 
^SL{Bd), Tsl{B^) and TsLiM), to check the pattern of power corrections in semileptonic 
decays. 



3 Electromagnetic Penguins and Rare B Decays in 
the Standard Model 

The SM does not admit FCNC transitions in the Born approximation, which is obvious 
from the Lagrangian given at the very outset of these lectures. However, they are induced 
through the exchange of bosons in loop diagrams. We shall discuss representative 
examples from several such transitions involving B decays, starting with the decay B — > 



Xs + 7, which has been measured by CLEO [Q. This was preceded by the measurement 
of the exclusive decay B — > K* + 7 pl| : 



B{B^X, + -f) = (2.32 ± 0.57 ± 0.35) X 10-\ (113) 
B{B^K* + -f) = (4.5 ± 1.0 ±0.6) X 10^^ , (114) 
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yielding an exclusive-to-inclusive ratio: 

T(B ^ K* + -y) 

These decay rates test the SM and the models for decay form factors and we shall study 
them quantitatively. 

The leading contribution to 6 ^ s + 7 arises at one-loop from the so-called penguin 
diagrams and the matrix element in the lowest order can be written as: 

M{b^s +7) = ^^Y.V^bV*F2{^,)q^'e•'sa,,{mhR + m,L)h , (116) 
V2 ^ 

where Xi = m^/m^, and are, respectively, the photon four-momentum and po- 
larization vector, the sum is over the quarks, u, c, and t, and Vij are the CKM matrix 
elements. The (modified) Inami-Lim function F2{xi) derived from the (1-loop) penguin 
diagrams is given by p2|: 



Fo(x) 



X 



24(x - 1)4 



6a;(3x-2)loga;-(x-l)(8a;^ + 5a;-7) , (117) 



where in writing the expression for F2{xi) above we have left out a constant from the 
function derived by Inami and Lim, since on using the unitarity constraint these sum to 
zero. It is instructive to write the unitarity constraint for the decays 5 -|- 7 in full: 

Vt,vz + + KfeKTs = . (118) 

Now, since the last term in this sum is completely negligible compared to the others (by 
direct experimental measurements), one could set it to zero enabling us to express the 
one-loop electromagnetic penguin amplitude as follows: 

G 

M{b^s + 7) = ^ At {F2{xt) - F2{xc)) q^e''sa^,{mbR + msL)b . (119) 

The GIM mechanism is manifest in this amplitude and the CKM-matrix element 
dependence is factorized in At = VtbV^*. The measurement of the branching ratio for 
B ^ Xg + 'y can then be readily interpreted in terms of the CKM-matrix element product 
At/|Vcfe| or equivalently |Vts|/|Vcfe|. In the approximation we are using (i.e., setting A„ = 
0), this is equivalent to measuring \Vcs\- For a quantitative determination of iV^sl/lKftl, 
however, QCD radiative corrections have to be computed and the contribution of the 
so-called long-distance effects estimated. We proceed to discuss them below. 



3.1 The effective Hamiltonian for B Xg-j 

The appropriate framework to incorporate QCD corrections is that of an effective theory 
obtained by integrating out the heavy degrees of freedom, which in the present context are 
the top quark and bosons. This effective theory is an expansion in 1/m^ and involves 
a tower of increasing higher dimensional operators built from the quark fields [u, d, s, c, b), 
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photon, gluons and leptons. The presence of the top quark and of the bosons is 
reflected through the effective coefficients of these operators which become functions of 
their masses. The operator basis depends on the underlying theory and in these lectures 
we shall concentrate on the standard model. The basis that we shall use is restricted to 
dimension-6 operators and it closes under QCD renormalization. The operators which 
vanish on using the equations of motion are not included. The effective Hamiltonian "He// 
given below covers not only the decay 6 — > s + 7, in which we are principally interested 
in this section, but also other processes such as b ^ s + g and b ^ s + qq. 

It is to be expected in general that due to QCD corrections, which induce operator- 
mixing, additional contributions with different CKM pre-factors have to be included in 
the amplitudes. Thus, QCD effects alter the CKM-matrix element dependence of the 
decay rates for both B ^ Xg + 'y and (more importantly) B — >■ X^ + 'J- However, with the 
help of the unitarity condition given above, the CKM matrix dependence in the effective 
Hamiltonian incorporating the QCD corrections for the decays i? — > + 7 factorizes, 
and one can write this Hamiltonian as |^: 

AC ^ 

n,ff{b - s + 7) = — ^v;:\4,5:a(/i)o.(/x), (120) 

where the operator basis is chosen to be (here /i and u are Lorentz indices and a and jS 
are colour indices) 

01 = {sLalt.bLa){cLpYcLf3), (121) 

02 = isLal^.bLf3)icL(37''CLa), (122) 

03 = (sLalf^bLa) E i^lL^l^ ^iL/s) , (123) 

q=u,d,s,c,b 

= (sLalM E (QL^rqia), (124) 
q=u,d,s,c,b 

O5 = {sLalt^bLa) E imi^m)^ (125) 
q=u,d,s,c,b 

Oe = (sLalt^bL/s) E iQRf^l^(lRa), (126) 

q=u,d,s,c,b 

= T^MsLaa^ubRa)F'''', (127) 
^7 = T^rn,{sR^a^,bLa)F'''', (128) 



^8 = Tt^MsLaT:f,a,,bRp)G''>^', (129) 

^8 = Tl^rn,{sR^T^^a,jLp)G^^\ (130) 

where e and gs are the electromagnetic and the strong coupling constants, and F^u and 
G^^ denote the electromagnetic and the gluonic field strength tensors, respectively. We 

^Note that in addition to the penguins with the u-quark intermediate state there are also non- 
factorizing contributions due to the operators (uLaj'^bLa)isLi3j^iULp), which hke the w-quark contribution 
to the 1-loop electromagnetic penguins are proportional to the CKM-factor A„ = KisVy*^, and hence are 
consistently set to zero. 
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call attention to the explicit mass factors in OjIOj) and Os{Og), which will undergo 
renormalization just as the Wilson coefficients. The dominant contributions in the radia- 
tive decays B ^ Xg + 'j arise from the operators O2, and Og, whereas the operators 
Os, get coefficients through operator mixing only, which numerically are negligible. 



Historically, the anomalous dimension matrix was calculated in a truncated basis 
and this basis is still often used for the sake of ease in calculating the real and virtual 
corrections, though as we discuss below, now the complete anomalous dimension matrix 
is available |^ . 



The perturbative QCD corrections to the decay rate r(i? — > + 7) have two distinct 
contributions: 

• Corrections to the Wilson coefficients calculated with the help of the renor- 
mahzation group equation, whose solution requires the knowledge of the anomalous 
dimension matrix in a given order in a^. 

• Corrections to the matrix elements of the operators Oi entering through the effective 
Hamiltonian at the scale ^ = 0{miy). 

The anomalous dimension matrix is needed in order to use the renormalization group and 
sum up large logarithms, i.e., terms like a^{mw) log™'(mf,/M), where M = rrit or mw and 
m < n (with n = 0,1,2, ...). At present only the leading logarithmic corrections (m = n) 
have been calculated systematically in the complete basis given above |^ . 



Next-to-leading order corrections to the matrix elements are now available completely. 
They are of two kinds: 

• QCD Bremsstrahlung corrections b ^ s'-f + g, which are needed both to cancel the 
infrared divergences in the decay rate for B — + 7 and in obtaining a non-trivial 
QCD contribution to the photon energy spectrum in the inclusive decay B Xg+'j. 

• Next-to-leading order virtual corrections to the matrix elements in the decay b 
5-^7. 



The Bremsstrahlung corrections were calculated in 1 95] - [9? in the truncated basis and 



last year also in the complete operator basis [Q, which have been checked recently in 



98|P]. The higher order matching conditions, i.e., Ci^rnw), are known up to the desired 
accuracy, i.e., up to 0{as{Mw) terms [^]. Very recently, the next-to-leading order virtual 
corrections have also been calculated ||100|| . What still remains to be done is the calculation 
of 7''^'*, the next-to- leading order anomalous dimension matrix, which is the hardest part 
of this computation. This is clearly needed to get theoretical estimates of the branching 
ratio B{B — >■ X<. + 7) in the complete next-to-leading order. We discuss the presently 
available pieces to estimate B{B — > + 7). 

We recall that the Wilson coefficients obey the renormalization group equation 



'9 ^, ,d' 



C.(^^,9yU9n[j^,9y (131) 



"This paper as it stands has some errors which we expect to be rectified! 
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The QCD beta function (3{g) has been defined earher and ^{g) is the anomalous dimension 
matrix, which, to leading logarithmic accuracy, is given by 



lig) = 7o 



167r2 



(132) 



Here 70 is a 8 x 8 matrix given by |]9^, |101 
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where / is the number of active fiavours, and u (d) is the number of up (down) flavours 
in the effective theory. For the case at hand we have u = 2, d = 3 and f = u + d = 5. 
This matrix is given in the NDR dimensional regularization scheme which is also used in 
the calculations of the matrix elements of the operators discussed below. The difference 
between the NDR and the HV schemes lies in the seventh and eighth columns and the 



corresponding matrix in the HV scheme can be seen in ||129 



The solution of the renormalization group flow is obtained as 

C7(/i) = T'^D{n,mw)TC{mw) 



(133) 



where T diagonalizes the anomalous dimension matrix and D is a diagonal matrix, which 
is given for the relevant case by 



D = Diag[r/°■l^5^ r/-°•899^ r]'^/'', 7^^^'^\ 7^-^^l^\ ry-O-^^so^ ^o.4086^ ^6/23] 



where 



w 



The transformation of the anomalous dimension matrix into a diagonal form can then 
be performed. Assuming that the coefficients C3 ■ ■ ■ Ce vanish at the matching scale 
H = Mw, the result for the coefficients at the scale 11 reads 



(134) 
(135) 



Cifi) = C2(Mh/) -0.0112r/-0-«994 + ^^-12/23 _Q^14Q3^-0.4230 ^Q_QQ54^0.1456 ^^gg) 

6 
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0.4086 



CM 
CM 

CM 



- 0.071^/2=^ + 0.0509r/' 

C2{Mw) ( 0.0156r/"°-«99^ - %-^2/23 ^ o.l214r7~°-423o ^ o.0026r7°-i456 (137) 
V 6 



- 0.071477^/23 ^ o.098477°-4°^6 
C2iMw) (-0.0025r7-°-«99^ + 0.0117r7-°-^230 ^ o.0304r7°-i^56 _ o.0397r7°-^°«6 



:i38l 



C2(Mvi/) (-0.0462r7-°-^994 + 0.0239r/~°-4230 _ o.0112r/°-i45^ + 0.0335r/ 



0.4086 



(139) 
(140) 



Csifx) 



C,{M^W^''' + C8(Mh.)^ (r7^^/23 - ^^6/23 
+C2{Mw) (-0.0185r/-o-«99^ - 0.0714r/-i2/23 _ o.0380r/-°-4230 _ q^qq^^^^-i^^^ 

-0.4286r7^/23 - 0.64947/°-^°^^ + 2.29967/^^/^3 - 1.08807/^^/^3 
C^8(Mh.)V^/23 

+C2(Mm/) (-0.05717/-°-^^^^ + 0.0873r/-°-^230 ^ 0.0209?/' 
-0.91357/°-^°«^ + 0.8623r/i4/23) . 



(141) 



,0.1456 



The non-zero initial conditions in the SM are given at the scale M\y and read |9 



C2 = 


1 


CriMw) = 


1 

X 

2 


Cs{Mw) = 


1 

X 

2 



2x73 + 5x/12 - 7/12 3xV2 - x 



(x-iy 



{x - ly 



\nx 



xV4 - 5x/4 - 1/2 3x/2 



{x-iy 



+ 



X - ly 



\nx 



(142) 
(143) 

(144) 



and X = m^/M^. The numerical values for the Wilson coefficients at the scale fi = M\y 
("Matching Conditions") and at three other scales /i = 10.0 GeV, 5.0 GeV and 10.0 GeV 
are given in Table ^ Also, for subsequent discussion it is useful to define two effective 
Wilson coefficients Cf (fx) and Cf (fx) jTOTl : 



= 07 — ^ — Oe 
Cf = Cs + C, . 



(145) 



Their values are also given in Table 1. 



3.2 Real and virtual 0{as) corrections for the matrix element 

Now, we discuss the real and virtual 0{as) corrections to the matrix element for 6^5 + 7 
at the scale fi ~ rrib, which by themselves form a well-defined gauge invariant albeit 
scheme-dependent set of corrections. This scheme dependence will be cancelled against 
the one in the anomalous dimension 7^^-*, as discussed in the context of the dominant 
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y(i — Tfl\\r 


jJ, — iU.U Kjii V 


jjj — O.U Kjii V 


jJj — Z.O Kjii V 


Ci 


0.0 


-0.158 


-0.235 


-0.338 


C2 


1.0 


1.063 


1.100 


1.156 




0.0 


0.007 


0.011 


0.016 


04 


U.U 


— U.Ui ( 


— U.Uz4 


— U.Uo4 




0.0 


0.005 


0.007 


0.009 


Ce 


0.0 


-0.019 


-0.029 


-0.044 


C7 


-0.193 


-0.290 


-0.333 


-0.388 


Cs 


-0.096 


-0.138 


-0.153 


-0.171 


Cf 


-0.193 


-0.273 


-0.306 


-0.347 


r^ff 


-0.096 


-0.132 


-0.146 


-0.162 



Table 1: Wilson coefficients Ci{fi) at the scale yU = mw = 80.33 GeV ("matching condi- 
tions") and at three other scales, /x = 10.0 GeV, = 5.0 GeV and /i = 2.5 GeV, evaluated 
with two-loop /3-function and the leading-order anomalous- dimension matrix. The entries 
correspond to the top quark mass ?T77(mf°'' = 170 GeV (equivalently, mf'" = 180 GeV) 
and the QCD coupling constant as{rn\) = 0.117, both in the MS scheme. 



contributions to the non-leptonic decays of the B hadron earlier. The results presented 
here correspond to the NDR scheme. 

Recapitulating the essential steps, we recall that the Bremsstrahlung corrections in 
b ^ sj + g, calculated in - and [Q, were aimed at getting a non-trivial photon 
energy spectrum at the partonic level. In these papers, the virtual corrections to 6 ^ 57 in 
0{as) were included only by taking into account those virtual diagrams which are needed 
to cancel the infrared singularities (and also the coUinear ones in the limit rris 0) 
generated by the Bremsstrahlung diagram. The emphasis was on deriving the photon 
energy spectrum in ^ -|- 7 away from the end-point —>■ 1 and the Sudakov- 
improved photon energy spectrum in the region — > 1. Clearly, the left-out virtual 
diagrams shown in Fig. ^ do not contribute either to the Sudakov spectrum or to the 
region 7^ 1 at the parton level. They, however, do contribute to the overall decay rate 
in B ^ Xs + 'j. Recently, these additional virtual correction have been evaluated in ||100|| , 
neglecting the contributions of the small operators O^-Oq. This additional contribution 
reduces substantially the scale dependence of the leading order (or partial next-to-leading 
order) decay width T{B + 7), which previously was found to be substantial and 

constituted a good fraction of the theoretical uncertainty in the inclusive decay rate 



102, 101, 103, E 



We shall follow closely the derivation of the virtual corrections given in [|100|| . Concen- 
trating on the dominant operators O2, 0-j and Og, the contribution of the next-to- leading 
order correction to the matrix element part in 6 — > s + 7 can be expressed as follows: 

M=M2 + Mj + Ms (146) 

and the various terms (including appropriate counterterm contributions) can be summa- 
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Figure 3: Non- vanishing two-loop diagrams associated with the operator O2 calculated in 



100 1 . The fermions (6, s and c quark) are represented by solid lines. The wavy (dashed) 



line represents the photon (gluon). 



rized as ||100|| : 



with 



M2={sj\0,\b),ree^ U log ^ + 



416 



(147) 
(148) 



243 
+ 

+ 

+ 



{-833 + lUn^z^/^ 

1728 - ISOn^ - 1296C(3) + (1296 - 3247r2)L + 108L^ + 36L' 
648 + 727r2 + (432 - 2167r2)L + 36^^ 
"-54 - 847r2 + 1092L - 756^^] } 



(149) 



167r r ^ 



45 - 37r2 + 9L + 9L'^ 



z + 



37r2 + 9^2] z'^ + [28 - 12L] z^ }-150) 



Here, 3?r2 and denote the real and the imaginary part of r2, respectively, z = {mc/mi,Y 
and L = log(z). 
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The real and virtual corrections associated with the operator O7, calculated in 



96| , |60|| can be combined into a modified matrix element for 6 — > 57, in such a way that 

its square reproduces the result derived in these papers. This modified matrix element 
_^rnod ^g^^g pq] : 



with 



= {Sl\0,\b),ree ( ^ + ^ ( log y + 1 1 (151) 



?7=3 , rj = -iA-n') . (152) 



Finally, the result for A^g is ||100|: 



Ms = {S^\0r\b),ree ^ (^4 log ^ + ^8 j , (153) 

with 

4 = -Y , r8 = -^ (-33 + 27r2-6m) . (154) 

With the results given above, one can write down the amplitude M{b ^ S7) by 
summing the various contributions already mentioned. Since the relevant scale for a b 
quark decay is expected to be /i ~ rub, the matrix elements of the operators may be 
expanded around /i = up to order 0(a<j) and the next-to- leading order result can be 
written as: 



M{b^ S7) = D {s^\Or{mb)\b)tree (155) 



with 



D = Cfi,) + ^ {cr"i,hP'^^ log ^ + Cr"r-^ , (156) 

where the quantities '^^^^^^^ = £j + 8 6i^ are the entries of the (effective) leading order 
anomalous dimension matrix and the quantities ii and are given for i = 2, 7, 8 in eqs. 
( |148|JT49|) , (HI) and (|15|), respectively. The first term, Cf{ii), on the r.h.s. of Eq. ([T56|) 



has to be taken up to next-to-leading logarithmic precision in order to get the full next- 
to- leading logarithmic result, whereas it is sufficient to use the leading logarithmic values 
of the other Wilson coefficients in Eq. ( |156| ). 

As pointed out by Buras et al. ||101|| , the explicit logarithms of the form ag {rrib) log(mfe / yu) 
in Eq. ( |156| ) are cancelled by the /^-dependence of C^^^^^ (fi). Therefore, the scale depen- 
dence is significantly reduced by including the virtual corrections completely to this order. 

The decay width F^*^* which follows from Ai{b ^ 57) in Eq. ( p.55| ) reads 



virt _ ^h,pole ^F^t 



em 



327r^ 



F\D\^ , (157) 



where the terms of 0(0;^) in \D\ have been discarded. The factor F in Eq. ( |1571 ) is 

/ \ 2 

fnibifi = nib)] 8 as{mb) 



mb,pole J 3 TT 



(158) 
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and its origin lies in the explicit presence of irib in the operator O7. To get the inclusive 
decay width for b — >■ 57(5'), also the Bremsstrahlung corrections (except the part already 
absorbed above) must be added. The contribution of the operators O2 and Oj was 



calculated already in [ p5[] . 



3.3 Estimating long-distance effects in B ^ Xg + ^ 

In order to get the complete amplitude for B Xg + j one has to include also the 
effects of the long-distance contributions, which arise from the matrix elements of the 
four-quark operators in Tie//, {Xs'y\Oi\B) . We shall discuss such contributions in the 
exclusive radiative decays -B — ^ (p, cu) + 7 subsequently, which have been calculated in 
a more robust theoretical framework using QCD sum rules. This framework (likewise 
lattice-QCD) is not applicable in inclusive decays in a straightforward manner. Lacking 
anything better, we discuss phenomenological models used in the literature in estimating 
the long-distance amplitude in radiative decays B ^ Xg + •y. 

In calculating the matrix elements of four-Fermion operators, it is usual to invoke 
the hypothesis of factorization. In the present context, factorization is combined with the 
additional assumption of vector meson dominance, involving the decays B J2i Vi+Xg — 
7 + Xg, where Vi = J/ipjip' , ... | ]104[ - | ]108|| . It should be remarked that non-leptonic 



decays, such as B ^ {J/ijj,ip' + Xg), by themselves are not under the quantitative control 
of the factorization-based framework due to the presence of significant non-factorizing 
pieces in the non-leptonic amplitudes - perhaps a sign of the breakdown of local quark- 
hadron duality (?). Phenomenologically, however, data on two-body and quasi two-body 
decays are well accounted for in terms of the Bauer-Stech-Wirbel effective parameters 
ai = Ci/Nc + C2 and 02 = Ci + C2/Nc [tL09|| . The details of this analysis and a host of other 
related decay modes can be seen in | |1 1U| | . Concerning the use of vector meson dominance. 



one has to ensure that the resulting amplitude A4{B ^ Xg + 'f) remains manifestly gauge 
invariant. In the present model this amounts to discarding the longitudinal polarization 
contribution in the non-leptonic decays B {J/tpjtp' , ...) + Xg, which in fact dominates 
the decay widths | |11(J| |, and keeping only the smaller contribution from the transverse 
polarization of J/ip, .... Following ||106| , |107|| , one can write the decay amplitude as: 



M{b^ sJ/ij)T = %a2V,,V: f"^^'^'"'^^ -sa^%l + i,)hq,el{q) , (159) 

where is defined as {'>P{(l)\c'')^c\Q) = —ig^{^(f)€-'^^{q). For the decays under consideration 
one needs the value of 02, which has been found to be \a2\ = 0.24 ± 0.04 | |110|| . One also 
needs to evaluate the coupling constant gvi.'f') at the point = 0. From leptonic decays 
of vector mesons, one gets, however, gviq^ = My). It has been remarked in literature, 
in particular in [|106| , p.07|| , that using gv{q^ = 0) = gv{,q^ = My) would substantially 



overestimate the long-distance contribution due to the expected dynamical suppression 
of the effective coupling gviq"^), as one extrapolates to the point = 0. Taking all this 
into account, an estimate of the long-distance amplitude from the intermediate J/ip state 
is |T06i : 

MLD{b sJ/i, 57) = p^a2VtbV:g{\e^^^P^)-sa^^{l + ^,)hF,, , (160) 

Zy Z o m ji^mb 
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where use has been made of the ■i/' to 7 conversion vertex, 

< 0\eJ^J^p >= (2e/3) < 0\c-f^c\i: >= -(2e/3)z(7^(0)e^ , (161) 

and we must use the value gj/ip{0) for the conversion couphng constant. Including all the 
(cc) resonances and the short distance contribution J^sd, the two-body decay amplitude 
(6 — s'-f) can be written as 

M{b^ .7) = -^^V,,V;A^m,Dif^) - a,\ ^ ^^^]-sa^^{l + ^,)hF,, , (162) 

where ipi represents the following vector cc resonant states: ipi'^S), ■ip{2S), '?/'(3770), 
■?/'(4040), ip{4160), and ■?/'(4415), and D is the function given earlier. Since in the LD- 
amplitude, gluon Bremsstrahlung is neglected, it is consistent to neglect this also in the 
SD-piece. Taking this estimate as giving the right order of magnitude for the long-distance 
contribution, Deshpande et al. ||106|| conclude that such long-distance effects can be as 



large as 10%. Other estimates, in particular by Golowich and Pakvasa ||107|| , lead to 
an even smaller long-distance contribution. Clearly, one can not argue very conclusively 
if such estimates are completely quantitative due to the assumptions involved. In fu- 
ture, one could improve these estimates by using data from HERA on elastic J/ip—, and 
■j/^'-photoproduction to get gj/^ and (?,/,' (0) directly, reducing at least the extrapolation 
uncertainties involved in the presently adopted procedure of extracting these coupling 
constants from the leptonic decay widths of each state and extrapolating to the point 
= using an Ansatz. 

3.4 Estimates of B{B + 7) in the Standard Model 

In the quantitative estimates of the SM branching ratio B{B — > + 7) given below 
we have neglected the LD-contributions. This is an assumption, which as we discussed 
above, we do not expect to work much better than 0(10%). It is theoretically preferable 
to calculate this quantity in terms of the semileptonic decay branching ratio 

B{B ^ X,7) = fJ^^l±lA^th s{B ^ Xlv,) , (163) 

J- SL 

where, the leading-order QCD corrected Tsl has been given earlier. The leading order 
power corrections in the heavy quark expansion, discussed in the context of the semilep- 
tonic decay rate, are identical in the inclusive decay rates for B — > + 7 and B Xiui, 
entering in the numerator and denominator in the square bracket, respectively, and hence 
drop out H, |0[. 



The error on the branching ratio B{B ^^7) comes from four different sources. We 
list them below: 

1. Amt'. The present value of is = 175 ± 9 GeV [Q, which is usually interpreted 
as the pole mass although this is not unambiguous. With this the running top 
quark mass in the MS scheme is rnl = 166 ± 9 GeV. However, in Eq. (|156|) , there 
is no distinction made between mf"'*^ and mj. To take both the theoretical and 
experimental errors into account, we take rrit = 170 ± 15 GeV. This leads to an 
error of ±5% in B{B —>■ Xg'j). 



39 



2. Ayu: The scale dependence is now reduced thanks to the work done in [ |10U| ]. It is 
usual to estimate the residual /i-dependence by taking a range m6/2 < < 2mb. 
This is plotted in Fig. ^ (solid curves), which shows that it yields an error of ±6% 
on B{B ^Xs + -f). 

3. Errors from the extrinsic parameters, 

(i) A(mb) (equivalently from A{mc/mb)), 

(ii) A{as{mz) (equivalently AA5, the uncertainty on the QCD-scale parameter), 

(iii) A{BRsl), the experimental uncertainty on the semileptonic branching ratio. 
Taking rric/mb = 0.29±0.02, B{B Xiu^) = (10.4±0.4)%, a,(mz) = 0.117±0.006 
[|], an uncertainty of ±12% was estimated on B{B ^ + 7) in [jBO . 



4. Last, but by no means least, there is at present an incalculable error due to the fact 
that Cy-^ is not yet available in the next-to-leading order. 

In view of this missing piece, it is not possible to give a completely corrected NLL predic- 
tion for B{B —>■ Xg + '-/) at present. In what follows, we shall replace C^^ by its leading 
log value, which yields the branching ratio 



B{B ^ X, + 7) = (3.20 ± 0.30 ± 0.38) x 10"^ (164) 

where the first error comes from the combined error on Am^ and A//, as can be seen in 
Fig. H, and the second from the extrinsic source. This estimate is subject to an additive 
renormalization due to the missing NLL anomalous dimension piece and the somewhat 
less tractable contribution of the LD-effect, discussed in the previous subsection but not 
included in the numerical estimates of the branching ratio (|164|) . However, if the NLL 



anomalous dimensions of the four-quark operators Oi, ...Oq, which are known are any 
useful guide for estimating Cj-^ , we do not expect a substantial change of this branching 
ratio in the complete NLL. Of course, the big challenge is to show that this indeed is the 
case. In view of this, we propose to add an additional theoretical uncertainty of order 
±10% on B{B ^ X, ± 7), yielding B{B ^ X, ± 7) = (3.20 ± 0.30 ± 0.38 ± 0.32) x 10"^ 
Combining the theoretical errors in quadrature gives 

B{B ^ X, ± 7) = (3.20 ± 0.58) x 10~^ (165) 

which we consider is a fairly robust estimate of the SM branching ratio at present. This can 
be improved only after the two estimated pieces are replaced by the actual calculations. 
The SM branching ratio B{B — > X^ ± 7) is compatible with the present measurement 
B{B ^ X, ± 7) = (2.32 ± 0.67) x 10"^ Q. Expressed in terms of the CKM matrix 
element ratio, one gets 



Wt 



ts 



IK 



0.85±0.12(expt) ±0.10(th), (166) 



cb\ 



which is within errors consistent with unity, as expected from the unitarity of the CKM 
matrix. 
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Branching ratio for b s^{g) calculated in [p.00|| with the parameters 



Figure 4: 

\Vts\/\V,t\ = 1, \Vtt\ = 1 
are explained in the text. 



m 



pole 



4.8 GeV and mc/mi, = 0.29. The different curves 



3.5 Photon energy spectrum in B ^ Xg + ^ 

The two-body partonic process b s'y yields a photon energy spectrum l/{r)dr{b —>■ 
57) = 6{1 — x), where the scaled photon energy x is defined as = {ml — ml) / {2 nib) x . 
Perturbative QCD corrections, such as 6 — *• 57 + (7, give a characteristic Bremsstrahlung 



,or X 



spectrum in x in the interval [0, 1] peaking near the end-points, 
and E'^ — s> (or X — s> 0), arising from the soft-gluon and soft-photon configurations, 
respectively. Near the end-points, one has to improve the spectrum obtained in fixed 
order perturbation theory. This is usually done in the region x — 1 by isolating and 
exponentiating the leading behaviour in ctem^s (/^)™'log^"(l — x) with m < n, where fx is 
a typical momentum in the decay B ^ Xg + The running of is a non-leading effect, 
but as it is characteristic of QCD it modifies the Sudakov-improved end-point photon 

[El- 



energy spectrum 



compared to its analogue in QED 



As long as the s- 

quark mass is non-zero, there is no coUinear singularity in the spectrum. However, parts 
of the spectrum have large logarithms of the form log(m^/m^), which are important 
near the end-point x —>■ but their influence persists also in the intermediate photon 
energy region and they have to be resummed |]60| , |115 . 



It has been observed in a number of papers [|112| , |116| , rT7 |, that the x-moments (scaled 
photon energy) in B ^ Xg + 'y and those involving lepton energy in the decay B X„£z/£ 
are related. Defining the moments as: 



yuW(S^X, + 7) 
M^''\B^XJue) 



r Jo dx 

rfxx"-^(-^^^ 
/o dx Te dx ' 



(167) 
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Jo dx 

The moments A^*^") have been worked out in the leading non-trivial order in perturbation 
theory and the results can be expressed as: 

_ 1 + ^C^(Alog2n + Slogn + const.) (168) 
where Cp = 4/3, the leading coefficient is universal with A = —1 ||114|| , and the non- 



leading coefficients are process dependent; B = 7/2 and B = 31/6 for B — > X^+'y 



and B — > Xuii'e, respectively. Measurements of the moments could eventually be used to 
relate the CKM matrix element Vts and Vub- 

How about calculating the photon energy spectrum in B ^ Xg + 'j and the lepton 
energy spectrum in i? ^ Xiue completely in QCD? Concentrating on B Xg + 'y, it 
is known that there is a gap between the end-points of the physical spectrum (^E^"-^ = 
{m\ — m\) / {2mB)), where X = {rriK + ^tt), and the partonic spectrum (^E^°'^ = {ml — 
m^)/(2mb)), and hence in the parton-model description there is a window in which 
remains empty. This gap can only be ffiled up by non-perturbative effects, which in their 
simple form can be attributed to the 6-quark motion in the B hadron. Phenomenological 
models were already proposed in the infancy of B physics which took into account non- 
perturbative smearing of the partonic spectra to ffil this gap. The question is if this effect 
can be computed in QCD proper. 

Attempts to calculate the photon and lepton energy spectra in the heavy quark expan- 
sion method lead to formal expressions which near the end-point are divergent [41|, [117 - 



118|| . The point is that near the end-point, the energy release for the light quark system for 
the final state is not of 0{mb) but of the order of the parameter A = MB—mb ~ 0{Aqqy)). 
Thus, the expansion parameter is no longer 1/ml and the operator product expansion - 
the backbone of the heavy quark expansion - breaks down. This divergent series in the 
effective theory has to be cleverly resummed and the distributions averaged over mo- 



mentum bins ||116|| to smooth the increasing number of derivatives of the delta function. 



5"(1 — x). This resummation allows to define an effective non-perturbative shape function 
|116|, |113|] , which can not be calculated in the effective theory, but one could use this con- 



cept advantageously to relate the energy spectra in the semileptonic decays B X^iui 
and B Xs+'y. It has been further argued in ||113|| that a Gaussian distribution, such as 



the one used in the Fermi motion model in [B^, 53] , with an appropriate definition of the 



effective momentum-dependent 6-quark mass, as prescribed by the heavy quark effective 
theory, is a good approximation to the shape function in 5 — >■ Xuii't and 5 — »• X + 7. 

We shall confine ourselves to the discussion of the photon energy spectrum, which 
combines the perturbatively computed spectrum, discussed in the previous section, with 
a model of the quark motion which fills the mentioned energy gap in the end-point region. 
We will then show a comparison of this spectrum with the measured photon energy 
spectrum in — >■ + 7 0|. In this model |21[], which admittedly is simplistic but has 



received some theoretical support in the HQET approach subsequently [^, |117|| , the b 



quark in B hadron is assumed to have a Gaussian distributed Fermi motion determined 
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by a non-perturbative parameter, pp, 



4 — 

^(p) = r- n exp( — ^) , p=\p\ (169) 

with the wave function normahzation /q°° dpp'^(f){p) = 1. The photon energy spectrum 
from the decay of the i?-meson at rest is then given by 



j'j"' dpp^ct>{p)j^{W,p,E,) , (170) 



dE^ Jo Lixv^ 

where Pmax is the maximally allowed value of p and is the photon energy spectrum 
from the decay of the 6-quark in flight, having a momentum-dependent mass W{p). This is 
calculated in perturbation theory taking into account the appropriate Sudakov behaviour 
in the end-point region at the partonic level. 

An analysis of the CLEO photon energy spectrum has been undertaken in |^| to 



determine the non-perturbative parameters of this model, namely mh{pole) and pp- The 
latter is related to the kinetic energy parameter Ai defined earlier in the HQET approach. 
The experimental errors are still large and the fits result in relatively small values; 
the minimum, xLin = 0.038, is obtained for pp = 450 MeV and mb{pole) = 4.77 GeV, 
in good agreement with theoretical estimates of the same, namely mb{pole) = 4.8 ± 0.15 
GeV H, 01 and p^ = -Ai/2 = 0.25 ± 0.05 GeV^ obtained from the QCD sum rules 



The central value, Ai = —0.4 GeV^, obtained from the photon energy spectrum in 
B Xs + 'J is also in good agreement with a recent determination of the same from 
the lepton energy spectrum in i? — ^ Xiui in which partially integrated spectrum and its 
first moment are analyzed in terms of the HQET parameters A and Ai using CLEO data, 
getting Ai = -0.35 ± 0.05 GeV^ and A = 0.55 ± 0.05 GeV 0. In Fig. § we have plotted 



the photon energy spectrum normalized to unit area in the interval between 1.95 GeV 
and 2.95 GeV for the parameters which correspond to the minimum (solid curve) and 
for another set of parameters that lies near the x^-boundary defined by = Xmin + 1- 



(dashed curve). Data from CLEO |90] are also shown. Further details of this analysis can 
be seen in Il60|. 



3.6 Inclusive radiative decays B + 7 

The theoretical interest in studying the (CKM-suppressed) inclusive radiative decays B —>■ 
X(i + 7 lies in the first place in the possibility of determining the CKM-Wolfenstein 
parameters p and rj. The relevant region in the decays B X^ + 7 is the end-point 
photon energy spectrum, which has to be measured requiring that the hadronic system 
Xfi recoiling against the photon does not contain strange hadrons to suppress the large- i?^ 
photons from the decay S — + 7. Assuming that this is feasible, one can determine 
from the ratio of the decay rates B{B + 'y)/B{B ^ X^ + 7) the CKM-Wolfenstein 

parameters. This measurement was proposed in where the photon energy spectra 
were also worked out. In spirit, such an analysis is very similar to the already undertaken 
for the inclusive semileptonic decays in which the end-point lepton energy spectrum is 
solved for the ratio |\4fe|/|V^b|. Now, the end-point spectra in B — >• X-t-7 will be analyzed 
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m^„p = 170 GeV 



450 MeV, 
310 MeV, m. 



300 MeV 




Figure 5: Comparison of the normalized photon energy distribution using the CLEO data 
90| corrected for detector effects and theoretical distributions from , both normalized 



to unit area in the photon energy interval between 1.95 GeV and 2.95 GeV. The solid 
curve corresponds to the values with the minimum x^, {■mq,pp) = {0, 450 MeV), and the 
dashed curve to the values (mg,p^) = (300 MeV, 310 MeV). 



in terms of the ratio |Vfd|/|Vts|. Of course, the experimental issues involved in the two 
cases are quite different and measuring B X^+j is lot more challenging than measuring 
B ^ XJue. 

In close analogy with the S — » + 7 case discussed earlier, the complete set of 
dimension-6 operators relevant for the processes b ^ d'j and b — >■ d'yg can be written as: 

Heffib ^d) = ^6 E^.(-")0,(/i), (171) 



where = Vjb V*^ for j = t, c, u. The operators O-,-, j = 1, 2, have implicit in them CKM 



factors. In the Wolfenstein parametrization |Tj], one can express these factors as : 

iu = A\\p-ir^), i, = -AX\ ^, = -^„-^,. (172) 

We note that all three CKM-angle-dependent quantities are of the same order of magni- 
tude, 0(A'^), where A = sin6'c — 0.22. This is an important difference as compared to the 
effective Hamiltonian 7ie//(& — > s) written earlier, in which case the effective Hamiltonian 
factorizes into an overall CKM factor = VtbV^l. For calculational ease, this difference 
can be implemented by defining the operators Oi and O2 entering in T-Ceff{b d) as 
follows 



01 = - — {cLf3l'^bLa){dLa7iMCLfi) " -r{uLI3'l^hLo,){dLc,'l fiULp) , 

02 = --f{cLaYbLa){dLplf,CLp) - -^iuLaYbLa){dLI3lfiULl3) , 



(173) 



and the rest of the operators {Oj; j = 3. ..8) are defined like their counterparts Oj in 
T-Ceff{b — *• s), with the obvious replacement s — *• ci. With this definition, the matching 
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conditions Cj{mw) and the solutions of the RG equations yielding Cj{n) become identical 
for the two operator bases Oj and Oj. 

It has been explicitly checked in the 0(as) calculations of the decay rate and photon 
energy spectrum involving b ^ dj and b dg'y transitions that the limit rUu — ^ 
for the decay rate T{B + 7) exists [^. This implies that the rates and E-y- 



spectrum are free of mass singularities. In a recent paper by Greub, Hurth and Wyler 



lOOj , the finiteness proof has been extended to the NLL order, in which the virtual 
corrections to the matrix elements are calculated completely. From this it follows that 
there are no logarithms of the type aemO^s log(m^/m^) ||108| , |100|| . Some papers, estimating 



LD-contributions in radiative B decays, seem to contradict this by assuming light-quark 
contributions which have such spurious log-dependence. There is no calculational basis for 
this assumption. On the other hand, as far as the dependence of the decay rate and spectra 
on the external light quark masses is concerned, one encounters logarithms of the type 
aemas[{^ + {'^-xf)/x] log{ml/m1) (for b sg-f) and aemas[(l + (1 log(m^/m^) 
(for b —>■ dg'y), which are important near the soft-photon (x —>■ 0) region and must 
also be exponentiated ||115 . 



The essential difference between T{B + 7) and T{B Xd + 7) lies in the 

matrix elements of the first two operators Oi and O2 (in Heff{b —>■ s)) and Oi and O2 (in 
Heff{b d)). The derivation of the inclusive decay rate and the final-state distributions 
in B —>■ Xd + 7 otherwise goes along very similar lines as for the decays B ^ Xs + j . 
The branching ratio B{B Xd + 7) in the SM can be written as: 



|2 



B{B^Xd + j)=D^\^t\ 

where the functions Di depend on the parameters mt^mi^^mc-, ^i-, as well as the others we 
discussed in the context of B{B ^ + 7) in the previous section. For the central values 
of these parameters , one gets Di = 0.21, D2 = 0.17, D3 = 0.03, = 0.10 m. This 



analysis has to be updated taking into account the complete NLL virtual calculations in 
|100|| , and so the numbers being quoted for Di are subject to some change. To get the 



inclusive branching ratio the CKM parameters p and rj have to be constrained from the 
unitarity fits. Taking the parameters from a recent fit, one gets 5.0 x 10"'^ < |^(| < 1.4 x 
10^^ (at 95% C.L.) ||120|| , yielding an order of magnitude uncertainty in B{B — > Xd + 7) 
- hence the interest in measuring this branching ratio. Taking the central values of the 
fitted CKM parameters discussed later in these lectures A = 0.8, A = 0.2205, 77 = 0.34 
and p = -0.07 one gets 

-B(5^Xd + 7) = (1.7±0.85) X 10-^ (175) 

which is approximately a factor 10 — 20 smaller than the CKM-allowed branching ratio 
B{B ^ X, + 7), measured by CLEO H. 



45 



3.7 Estimates of B{B 
parameters 



y + 7) and constraints on the CKM 



Exclusive radiative B decays B ^ V + •y, with V = K*,p,uj, are also potentially very 
interesting from the point of view of determining the CKM parameters ||121|| . The extrac- 
tion of these parameters would, however, involve a trustworthy estimate of the SD- and 
LD-contributions in the decay amplitudes. 

The SD-contribution in the exclusive decays Bd) — > {K*, p) + Bd ^ uj + •y and 
the corresponding Bs decays, Bg (0, K*) + 7, involve the magnetic moment operator 
Oj and the related one obtained by the obvious change s ^ d, O7. The transition form 
factors governing the radiative B decays B V + y can be generically defined as: 



{V,X\-^Pa,,q%\B) 



WrnP 



(176) 



Here \^ is a vector meson with the polarization vector e^'^-', V = p,uj, K* or (p; B is a generic 
5-meson B^, B^ or Bg, and ip stands for the field of a light u,d or s quark. The vectors 
Pb, Pv and q = Pb — Pv correspond to the 4- momenta of the initial i?-meson and the 
outgoing vector meson and photon, respectively. In ( |1 761 ) the QCD renormalization of the 
'^o'^i,q'^b operator is implied. Keeping only the SD-contribution leads to obvious relations 
among the exclusive decay rates, exemplified here by the decay rates for {Bu, Bd) p + 7 
and {Bu, Bd) ^ K* + r 



r{{B^,Bd) 


-^p + 7) 






T{{B^,Bd) - 


-^K* + y) 







u,d 



K. 



u,d 



IK 



td\ 



Wt 



ts\ 



(177) 



where is a phase-space factor which in all cases is close to 1 and 



^B,^K'y\ 



l2 



(178) 



is the ratio of the (SD) form factors squared. The transition form factors Fs are model 
dependent, and since the exclusive-to-inclusive ratio Rk* has been measured, one could 
use data to distinguish among models. This aside, the ratios of the form factors, i.e. Hi, 
should be more rehably calculable as they depend essentially only on the SU(3)-breaking 
effects. If the SD-amplitudes were the only contributions, the measurements of the CKM- 
suppressed radiative decays {Bu, Bd) P+7, Bd uj + y and Bg — > -ft'* +7 could be used 
in conjunction with the decays (i?„, Bd) K* + 7 to determine the CKM parameters. 
The present experimental upper limits on the CKM ratio |V^td|/|Vjs| from radiative B 



decays are indeed based on this assumption, yielding at 90% C.L. [^ : 



Vtd 

Vts 



< 0.64-0.76 



(179) 



depending on the models used for the SU{3) breaking effects in the form factors ||121|, \12'2 



The possibility of significant LD-contributions in radiative B decays from the light 
quark intermediate states has been raised in a number of papers [|104|] - ||108|| . Their 
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amplitudes necessarily involve other CKM matrix elements and hence the simple fac- 
torization of the decay rates in terms of the CKM factors involving \Vtd\ and \Vts\ no 
longer holds thereby invalidating the relationships given above. The discussion about the 
LD-contribution that we presented in the context of the inclusive decays B ^ Xg + 'y 
applies in its essence also for the exclusive decays, such as 5 — > K* + 7, with appropriate 
modifications. In line with this discussion, we shall assume that the LD-contributions are 
small also in exclusive decays B — > K* + 7. In what follows, we discuss some charged 
and neutral exclusive 5-decays, B^ — > p^7 and -B-decays B {p^,u)'j, involving the 
CKM-suppressed transitions and estimate the SD- and LD-contributions. 

The LD-contributions in B ^ V + 'y are induced by the matrix elements of the four- 
Fermion operators Oi and O2 (likewise Oi and O2). Estimates of these contributions 
require non-perturbative methods. This problem has been investigated recently in [|123| , 
|124|| using a technique [p.25|| which treats the photon emission from the light quarks in a 



theoretically consistent and model-independent way. This has been combined with the 
light-cone QCD sum rule approach to calculate both the SD and LD — parity conserving 
and parity violating — amplitudes in the decays Bu^d — > p{uj) + 7. To illustrate this, 
we concentrate on the 5^ decays, B^ ^ + 7 and take up the neutral B decays 
Bd — > p{uj) -|- 7 at the end. The LD-amplitude of the four-Fermion operators Oi, O2 is 
dominated by the contribution of the weak annihilation of valence quarks in the B meson. 
It is color-allowed for the decays of charged B^ mesons, as shown in Fig. ^ where also 
the tadpole diagram is shown, which, however, contributes only in the presence of gluonic 
corrections, and hence neglected. In the factorization approximation, one may write the 
dominant contribution in the operator O2 (here O2 is the part of O2 with the CKM factor 
in Eq. (^) 

(pimB) = {p\dr^u\o){-f\ur^^b\B) + {p^\dr^u\o){o\ur^b\B) , (iso) 

and make use of the definitions of the decay constants 

{0\uT^b\B) = tpjs, 

{p\dT,u\0) = e^fWp/,, (181) 

to reduce the problem at hand to the calculation of simpler form factors induced by vector 
and axial-vector currents. 

The factorization approximation assumed in | ]123| , |124| has not been tested experimen- 



tally in radiative B decays. From a theoretical point of view, non-factorizable contribu- 
tions belong to either the 0(0;^) (and higher order) radiative corrections or to contribu- 
tions of higher-twist operators to the sum rules. Their inclusion should not change the 
conclusions substantially. 

The LD-amplitude in the decay 5^ ^ + 7 can be written in terms of the form 
factors and F2, 

x{-i[g^''{q-p)-p^q^]-2Fi:{q^) + e^''"^p^qp-2F^{q^)]. (182) 
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Figure 6: Weak annihilation contributions in Bu P7 involving the operators 0[ and 
O2 denoted by with the photon emission from a) the loop containing the b quark, b) the 
loop containing the light quark, and c) the tadpole which contributes only with additional 
gluonic corrections. 



Again, one has to invoke a model to calculate the form factors. Estimates from the 
light-cone QCD sum rules give |124|: 



F[/Fs = 0.0125 ± 0.0010 , F^/Fg = 0.0155 ± 0.0010 , 



(183) 



where the errors correspond to the variation of the Borel parameter in the QCD sum 
rules. Including other possible uncertainties, one expects an accuracy of the ratios in 
( |183|) of order 20%. Estimates of Fs in the normalization of Eq. ( |176|) range between 
FsiB 



K*-f) = 0.31 (Narison in |T22[) to Fs{B K*-f) = 0.37 (Ball in ^^), with a 
typical error of ±15%, and hence are all consistent with each other. This, for example, 
gives Rk* = 0.16 ± 0.05, using the result from ||121|| , which is in good agreement with 
data. 

Returning to the discussion of the LD-contribution, we note that the parity-conserving 
and parity- violating amplitudes turn out to be numerically close to each other, Ff" ~ F2 = 
Fl, hence the ratio of the LD- and the SD- contributions reduces to a number |124 



hort 



R 



Using C2 = 1.10, Ci = 
^ = 5 GeV) H gives: 



-0.235, Cj^ = —0.306 from Table ^ (corresponding to the scale 



^L/S 



47r2mp(C2 + Ci/iVe) F^ 



Bu^pf 



rribCj 



-0.30 ±0.07 



:i85) 



which is not negligible. To get a ball-park estimate of the ratio Along/ -^short, we take the 
central values of the CKM matrix elements, = 0.9744 ±0.0010 H], \Vtd\ = (1.0 ±0.2) x 



10" 



\Vcb\ = 0.039 ± 0.004 and |14fe/K 



cb 



•Along/ Ashort] ^ 



0.08 ±0.02 g, yielding, 

\VubVud\ 



I p-««^P7| 



\VuVbt\ 



10% 



(186) 
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Thus, the CKM factors suppress the LD-contributions. 

The analogous LD-contributions to the neutral B decays — > pj and B^ — >■ U'j are 
expected to be much smaller, a point that has also been noted in the context of the VMD 
and quark model based estimates ||104|| . In the present approach, the corresponding form 
factors for the decays fi^ p^{uj)'~f are obtained from the ones for the decay B^ —* p^7 
discussed above by the replacement of the light quark charges e„ — > e^, which gives the 
factor —1/2; in addition, and more importantly, the LD-contribution to the neutral B 
decays is colour- suppressed, which reflects itself through the replacement of the factor ai 
by 02- This yields for the ratio 



-"-L/5 



-0.13 ±0.05, 



(187) 



where the numbers are based on using 02/01 = 0.27 ± 0.10 [ |110|| . This would then yield 



at most R 



L/S 



R 



L/S 



0.05, which in turn gives 



'^long 
'^short 



< 0.02. 



(188) 



Even if this underestimates the LD-contribution by a factor 2, due to the approximations 
made in [p.23|, |124|| , we conclude that it is quite safe to neglect the LD-contribution in the 



neutral 5-meson radiative decays. 

Restricting to the colour-allowed LD-contributions, the relations, which obtains ignor- 
ing such contributions (and isospin invariance). 



T{Bu ^ p+7) = 2 T{Bd ^ p°7) = 2 T{Bd ^ ^7) 



(189) 



get modified to 



r(i?„ ^ p7) T{Bu ^ p7) 



2r(5, ^ p7) 2T{B, 



1 + 2 ■ Rr/sV, 



L/S^ud- 



^7) 
-P) 



1 + rITs 



T/ T/* 
>P7 ^ub^ud 



VtbV;, 



p^ + r/^ 



ud I 



pY + -rf 



(190) 



where Rl/s = Rl/^^ ■ The ratio r(5„ rf)/2T{Bd -> P7)(= r(B„ ^ Pl)/2T{Bd ^ 
uj'-f)) is shown in Fig. |^ as a function of the parameter p, with 77 = 0.2, 0.3 and 0.4. This 
suggests that a measurement of this ratio would constrain the Wolfenstein parameters 
{p,r)), with the dependence on p more marked than on 77. In particular, a negative value 
of p leads to a constructive interference in B^ p7 decays, while large positive values of 
p give a destructive interference. This behaviour is in qualitative agreement with what 
has been also pointed out in ||105| . 

The ratio of the CKM-suppressed and CKM-allowed decay rates for charged B mesons 
gets modified due to the LD contributions. Following ||107|| , we ignore the LD-contributions 
in r{B — ^ K*'~f). The ratio of the decay rates in question can therefore be written as: 

r(5n ^ P7) 



r(5„ ^ K*j) 



P? + v'] 



1 + 2 ■ RL/sVud- 



p(l -p) -ri^ 
1 - p)2 + r/2 



+ iRL/s)X 



2 _ 

ud I 



p^ + 
1 - p)2 + r/2 



:i9i^ 
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Using the central value from the estimates of the ratio of the form factors squared 
Ku = 0.59 ±0.08 [ |121|| , we show the ratio (|191|) in Fig. || as a function of p for rj = 0.2, 0.3, 
and 0.4. It is seen that the dependence of this ratio is rather weak on r] but it depends on 
p rather sensitively. The effect of the LD-contributions is modest but not negligible, intro- 
ducing an uncertainty comparable to the ~ 15% uncertainty in the overall normalization 
due to the S'f/(3)-breaking effects in the quantity k^. 



\ 



Figure 7: Ratio of the neutral and charged i?-decay rates T{Bu — >■ p'y)/2T{Bd — > p7) 
as a function of the Wolfenstein parameter p, with r) = 0.2 (short-dashed curve), t] = 0.3 
(solid curve), and rj = 0.4 (long-dashed curve). (Figure taken from ||124|| .) 



Neutral 5-meson radiative decays are less-prone to the LD-effects, as argued above, 
and hence one expects that to a good approximation the ratio of the decay rates for 
neutral B meson obtained in the approximation of SD-dominance remains valid ||121||: 



T{Bd P7,^7) 
T{B ^ K*-f) 



n,\'[il-pf + r]' 



(192) 



where this relation holds for each of the two decay modes separately. It is a realistic hope 
that this relation is theoretically (almost) on the same footing in the standard model as 
the one for the ratio of the B^-B^ mixing-induced mass differences, which satisfies the 
relation 



120 



AM, 



Vt 



ts 



Vt 



td 



1 



A2[(l-p)2 + r/2 



(193) 



The hadronic uncertainty in this ratio is in the SU(3)-breaking factor Kgd = {fBs^BjfB^BsJ 
which involves the pseudoscalar couphng constants and the so-called bag constants. This 



50 



0.4 



Figure 8: Ratio of the CKM-suppressed and CKM-allowed radiative S-decay rates 
T{Bu — > frf)/r{B — i> K*'y) (with B = Bu or Bd) as a function of the Wolfenstein pa- 
rameter p, a) with 7] = 0.2 (short-dashed curve), rj = 0.3 (sohd curve), and t] = 0.4 
(long-dashed curve). (Figure taken from ||124] ].) 



quantity is estimated as Kgd = 1-35 ± 0.25 in the QCD sum rules and lattice QCD ap- 
proaches. (For details and references, see [|120|| ). The present upper limit for the mass- 
difference ratio AMs/AMd > 12.3 at 95 % C.L. from the ALEPH data [12§ provides 
a better constraint on the CKM parameters, yielding |Vj(i|/|\4s| < 0.35, than the corre- 
sponding constraints from the rare radiative decays B (p, u) + 7, which give an upper 
limit of 0.75 for the same CKM-ratio. We expect experimental sensitivity to increase 
in both measurements, reaching the level predicted for this ratio in the standard model, 
|^td|/|Vfs| = 0.24 ± 0.05 [|120|| , in the next several years in the ongoing experiments at 
CLEO, SLC and Tevatron, and the forthcoming ones at the B factories, HERA-B and 
LHC. 



Finally, combining the estimates for the LD- and SD-form factors in ||124|| and ||121 
respectively, and restricting the Wolfenstein parameters in the range —0.4 < p < 0.4 and 



0.2 < T] < 0.4, as suggested by the CKM-fits ||120|| , we give the following ranges for the 
absolute branching ratios: 



(194) 



B{B^^fn) = (1.9±1.6) X 10-^ 

B{Bd p-i) ~ B{Bd uj-f) = (0.85 ± 0.65) x lO"*' 



where we have used the experimental value for the branching ratio B{B — > K* + 7) pl| , 
adding the errors in quadrature. The large error reflects the poor knowledge of the CKM 
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matrix elements and hence experimental determination of these branching ratios will put 
rather stringent constraints on the Wolfenstein parameter p. 

Summarizing the effect of the LD-contributions in radiative B decays, we note that 
they are dominantly given by the annihilation diagrams. QCD sum-rule-based estimates 
are encouraging in that they strengthen the hope that such contributions are modest in 
exclusive radiative B decays, in particular in the neutral i?-decays B^ [p^^uj) + 7, 
with B^ — »• p^7 modified by O(20)% from its SD-rate. This should be checked in other 
theoretically sound frameworks. Of course, forthcoming data on charged and neutral B- 
meson decays will be able to determine the LD-contribution directly. Presently available 
data in general suggest that the contribution of annihilation diagrams in B decays is not 
significant, as seen through the near equality of the lifetimes for the S^, B^ and 5° 
mesons and the near equality of the observed -B^ and B^ radiative decay rates. In terms 
of the operator product expansion, they all involve the four-Fermi operators, sandwiched 
between various inclusive and exclusive states. The matrix elements are process-dependent 
and have to be estimated as well as possible. Some illustrative examples were given here. 



3.8 Inclusive rare decays B Xsi~^£ in the SM 

The decays B Xsi~^i~ , with i = e, p, r, provide a more sensitive search strategy 
for finding new physics in rare B decays than for example the decay B Xs'y , which 
constrains only the magnitude of Cj^ . This experimental constraint has triggered a lot 
of theoretical investigations on beyond-the-SM frameworks, as can be judged from the 
incomplete list of papers in ||127| , |128|| . The sign of Cj-^ , which depends on the underlying 



physics, is not determined by the measurement of B{B X,, + 7). This sign in our 
convention is negative in the SM (see Table 0). However, it is known (see for example 



129|| ) that in supersymmetric models, both the negative and positive signs are allowed as 
one scans over the allowed parameter space of this model. The determination of the sign 
of Cy'^ is an important matter as this will impose further constraints on the parameters 
of many models. It will also test the prediction of the SM, by the same token. 

Continuing this discussion, we recall that a part of the amplitude for B — > Xs£^£^ 
involving the coupling of the virtual photon to the charged lepton pair depends on 
the effective Wilson coefficient C^'^ encountered in the electromagnetic penguin decays 
— > X<i + 7. For low dilepton masses, the differential decay rate for B Xsi'^i~ 
is dominated by this contribution. However, as we shall see below, the B Xgi'^i' 
amplitude in the standard model has two additional terms, arising from the two FCNC 
four-Fermi operators, Q which are not constrained by the B ^ Xg + 'y data. Calling 



their coefficients Cg and Cio, it has been argued in ||129|| that the signs and magnitudes 



of all three coefficients Cj , Cg and Cio can, in principle, be determined from the decays 
B Xg + J and B — ^ Xsi~^£~ . The coefficient Cg, which governs the strength of the 
chromomagnetic moment transition, can, in principle, also be determined by measuring 
the photon energy spectrum in i? — > + 7 in low-to-intermediate photon energy region. 
The SM-based rates for the decay b — > si~^i~ , calculated in the free quark decay 



**This also holds for a large class of models such as MSSM and two-Higgs doublet models but not 
for all. In LR symmetric models, for example, there are additional FCNC four-Fermi operators involved 
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approximation, have been known in the LO approximation for some time [Q. The LO 
calculations have the unpleasant feature that the decay distributions and rates are scheme- 
dependent. The required NLO calculation is in the meanwhile available, which reduces the 
scheme-dependence of the LO effects in these decays [ |13CI| , |131|| . In addition, long-distance 
(LD) effects, which are expected to be very important in the decay B X^i^i^ | |132 |, 



have also been estimated from data on the assumption that they arise dominantly due to 
the charmonium resonances J/ip and ip' through the decay chains B —>■ XsJ/iIjIi/j' , ...) 
Xsi~^i~. The effect of these resonances persists even far away from the resonant masses 
deforming the short-distance based distributions appreciably [|132|| . Likewise, the leading 
{l/rrih^) power corrections to the partonic decay rate and the dilepton invariant mass 
distribution have been calculated with the help of the operator product expansion in the 
effective heavy quark theory ||133|| . A theoretically complete description including the 



improved perturbative treatment of the decay B — >■ Xsi~^i~ , the LD-effects and power 
corrections mentioned above, to the best of our knowledge, is still lacking. In what follows, 
we shall not take into account the power corrections but include the rest to get estimates 
for the decay rates and distributions in B ^ Xsi~^i~ . 

The amplitude for B —>■ Xsi~^i~ is calculated in the effective theory approach, which 
we have discussed earlier, by extending the operator basis of the effective Hamiltonian 
introduced in Eq. ( |120| ): 



Heffib ^s + r,b-^s + tr) = Heffib ^ s + 7) - -^VZVt, [C9O9 + CioOio] , (195) 



V2 



where the two additional operators are: 



a 



O9 = —Sal^PLbJl/, 

OL - 

Olo = -r'^.l^'PLbdl^.l^t . (196) 
47r 

The analytic expressions for C^irny^') and C\Q{rnw) can be seen in |129|| and will 



not be given here. We recall that the coefficient Cg in LO is scheme-dependent. However, 
this is compensated by an additional scheme-dependent part in the (one loop) matrix 
element of Og. We call the sum ^ which is scheme- independent and enters in the 
physical decay amplitude given below, 

M{b^s^tn = ^^v:y,,- 



(197) 



with 



Cf{s)^C,r^{s) + Y{s). (198) 
The function Y{s) is the one-loop matrix element of Og and is defined as: 

Y{s) = g{m„s){3Ci + C2 + 3C3 + C4 + 3C5 + Ce) 
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vis) 



--g{l,s) (4C3 + 4C4 + 3C5 + C6) 

-^giO,rs) (C3 + 3C4) 

+1{3C3 + C, + 3C5 + C,) 
y 

_^l(3Ci + C2-C3-3C4), 

1 + ^u^is) . 
n 



(199) 
(200) 



Here s = (p+ +p^)'^/ml is the scaled dilepton mass, with p± = {E±,p±) denoting four- 
momenta of i^, rhi are the scaled quark masses, rfii = mi/mh, and the function uj{s) 
represents the 0(0;^) correction from one-gluon exchange in the matrix element of Og, the 
analogue of which we have discussed in the context of the semileptonic B decays earlier, 
derived in [p5| : 



UJiS 



-2 4 2 5 + 4s 
-IT i^i2{s) msm(l — s) — - 



9 3 ^ ' 3 
2s(l + s)(l-2s 



3(1 - s)2(l + 2s) 



In s + 



3(1 + 2s) 
5 + 9s - 6s2 



ln(l - s) 



6(1 - s)(l + 2s) 



(201) 



The function g{z, s) includes the charm quark-antiquark pair contribution [ 130 , 131]: 
9{z, s) 



9 fi 



81 8 4 2,^ , f- : 



0(1 - y){\n ^ + _ ,7r) + 0(y - 1)2 arctan -yl= 



^(0,.) = ^-gln(-) 



4, . 4. 
- m s H — tiT 
9 9 



(202) 



(203) 



where y = Az / s. We recall from the discussion in [[130| , |131|| that ^ is dependent on the 
dimensional regularization scheme, with. 



(NDR) 
-1 (HV). 



(204) 



For numerical estimates, the Wilson coefficients Ci,...Cq and C^-^ are given in Table |] 
and we give here the value for Cg (in the NDR-scheme) and Cio for the central values of 
the parameters used in Table ^ (m^ = 170 GeV and Aqcd = 0.195 GeV): 



C^v/^«(^ = 5.0+1^ GeV) 



4 09+°-^^ 



1+5.0 
-2.5 

Cw{mw) = -4.32 . (205) 
Note that the coefficient Cio{mw) does not get renormalized by QCD corrections. 
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Following refs. | |129|| and |135| ], the differential decay rates in B Xsi~^i (ignoring 
lepton masses) are, 



dB{s) 

ds 



B 



a 



A? 



si 



Ms) 



(206) 



with u{s) 



's - (1 



'1 — ms)2], f{z) has been defined earlier as we discussed 



2a,(/i)/37r [(tt^ - 31/4)(1 - + 3/2], and 



ai[s,ms) 
02(5, rhs) 

^eff 



2\2 



-2r + s(l + mf) + (l-mf) 

-(1 + ml)s^ - (1 + 14m^ + m^)s + 2(1 + m^)(l - m^) 
(1 -m2)2 _ (1 + m2)s. 



2\2 



(207) 
(208) 
(209) 



Here '^{C^ ) represents the real part of C7 . A useful quantity is the differential FB 



asymmetry in the c.m.s. of the dilepton defined in refs. 135 | 



dA{s) 
ds 



dz 



^dB 

dz ' 



where z = cos 9, which can be expressed as: 



dAjs) 
ds 



-B. 



3a' 



47r2 /(m^ 



(210) 



(211) 



The Wilson coefficients Cj^^ , C^^^ and Cio appearing in the above equations can be deter- 
mined from data by solving the partial branching ratio S(As) and partial FB asymmetry 
^(As), where As defines an interval in the dilepton invariant mass ||129|| . 

There are other quantities which one can measure in the decays B Xsi^i^ to 
disentangle the underlying dynamics. We mention here the longitudinal polarization of 
the lepton in i? ^ Xsi~^i~ , in particular in i? ^ X^r+r^, proposed by Hewett [ I3ti| ]. In 



a recent paper, Kriiger and Sehgal ||137|| have stressed that complementary information 
is contained in the two orthogonal components of polarization [Pt, the component in 
the decay plane, and P/v, the component normal to the decay plane), both of which are 
proportional to m^/mt, and therefore significant for the r+r" channel. A third quantity, 
called energy asymmetry, proposed by Cho, Misiak and Wyler ||138||, defined as 



^ N{E,- >Ee+)- N{E,+ > E, 
N{Ee- > Ee+) + N{Ee+ > E^ 



(212) 



where N{Ei- > Ei+) denotes the number of lepton pairs where is more energetic 
than i~ in the 5-rest frame, is, however, not an independent measure, as it is directly 
proportional to the FB asymmetry discussed above. The relation is ||140|| : 



Ja{s) 



BxA. 



(213) 
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This is easy to notice if one writes the Mandelstam variable u(s) in the dilepton cm. and 
the S-hadron rest systems. 

Next, we discuss the effects of LD contributions in the processes B Xgi^i^ . Note 
that the LD contributions due to the vector mesons such as J/ip and as well as 
the continuum cc contribution already discussed, appear as an effective {sL'~^^hi){I'^^tj 
interaction term only, i.e. in the operator Og. This implies that the LD-contributions 
should change Cg effectively, C^ as discussed earlier is dominated by the SD-contribution, 
and Cio has no LD-contribution. In accordance with this, the function Y{s) is replaced 
by, 

Y{s) ^ Y'{s) = Y{s) + Fres(s), (214) 



where Iresl-s) is given as |135 



>^res(S) = — «: (3Ci + C2 + SCg + C4 + 3C5 + Cg) ^ ^.^ . ^ ^ 



(215) 

where k is a fudge factor, which appears due to the inadequacy of the factorization frame- 
work in describing data on i? ^ J /ipXg. Here we use k (3Ci + C2 + 86*3 + C4 + 3C5 + Cg) ■ 
+0.88 for the numerical calculation, which reproduces (in average) the measured branch- 
ing ratios for B — >■ J/tpX^ and B — >■ ip'Xg, after the contributions from the Xc states 



have been subtracted. This is consistent with the treatment of LD-effects in ||137|| and 
||139|| , where further discussions and theoretical uncertainties on the LD-contributions can 
be seen. The long-distance effects lead to significant interference effects in the dilepton 
invariant mass distribution and the FB asymmetry in i? ^ Xsi^i~ shown in Figs. ^ and 
p!0| , respectively. This can be used to test the SM, as the signs of the Wilson coefficients 
in general are model dependent. 

The leading (l/mg) power correction to the dilepton mass spectrum in i? ^ Xsi^i~ 
has been worked out in | ]141| | using heavy quark expansion. The correction was found to 
be positive and around 10% over a good part of the dilepton mass spectrum with the 
particular choice of the parameters Ai and A2 adopted in this paper. In the meanwhile, 
the estimate of the parameter Ai has changed (it has a negative value now) compared to 
what has been used in [|141j| . Apart from this, the analytic form of the power corrected 



result in [ 141 is also somewhat puzzling and remains to be verified as, on general grounds, 
one expects the heavy quark expansion to break down near the end-point of the dilepton 
mass spectrum, i.e., as s — s™"^, which is not suggested by the power-corrected spectrum 



reported in ||141|| . In the opinion of this author, the end-point spectrum requires integra- 
tion over a range of dilepton masses (or smearing) to be calculable in the heavy quark 
expansion method. A detailed discussion of this point and a new derivation of power 



corrections to the dilepton mass spectrum and FB asymmetry will be presented in ||140| 



Taking into account the spread in the values of the input parameters, /x. A, m^, and 
BsL discussed in the previous section in the context of B{B — + 7), we estimate the 
following branching ratios for the SD-piece only (i.e., from the intermediate top quark 
contribution only) ||140|| : 



B{B X,e+e-) = (8.4 ± 2.2) x 10"^ 
i3(5 ^ X,/i+/i") = (5.7 ± 1.2) X 10"^ 
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Figure 9: Dimuon invariant mass distribution in B ^ Xs/i+/i~ in the SM including next- 
to-leading order QCD correction. The dashed curve corresponds to the short- distance 
contribution only and the solid curve is the sum of the long-distance and short- distance 
contributions. (Figure taken from 



B{B ^ XsT+T-) = (2.6 ± 0.5) X 10" 



(216) 



where theoretical errors and the error on Bsl have been added in quadrature. These 
estimates are consistent with the results presented in [|131| , |137|| . The present experimental 
limit for the inclusive branching ratio in B ^ Xsi^i~ is actually still the one set by the 
UAl collaboration some time ago ||142|| , namely B{B Xsfi'^fJ^~) > 5.0 x 10^^. As far as 
we know, there are no interesting limits on the other two modes, involving X^e'^e" and 

XsT+T-. 

It is obvious from Fig. ^ that only in the dilepton mass region far away from the res- 
onances is there a hope of extracting the Wilson coefficients governing the short- distance 
physics. The region below the J/tp resonance is well suited for that purpose as the 
dilepton invariant mass distribution there is dominated by the SD-piece. Including the 
LD-contributions, following branching ratio has been estimated for the dilepton mass 



range 0.2 < s < 0.36 in 140 : 



B{B X,^V") = (1-3 ± 0.3) X 10" 



(217) 



with B{B XgC^e ) ~ B{B Xsfi^fi ). The FB-asymmetry is estimated to be in the 



range 10% - 27%, as can be seen in Fig. |T0|. These branching ratios and the FB asymmetry 
are expected to be measured within the next several years at BABAR, BELLE, CLEO, 
CDF, DO, and HERA-B. In the high invariant mass region, the short- distance contribution 
dominates. However, the rates are down by roughly an order of magnitude compared to 
the region below the J/ip-mass. Estimates of the branching ratios are of 0(10^''), which 
should be accessible at the LHC. 



57 



->Xs mu+ mu- 



0.4 - 



0.2 




-0.2 



-0.4 - 



Figure 10: FB asymmetry for B X^fJ^^fi^ in the SM as a function of the dimuon 
invariant mass including the next-to-leading order QCD correction. The dashed curve 
corresponds to the short-distance contribution only and the solid curve is the sum of the 
long-distance and short-distance contributions. (Figure taken from ||140|| .) 



while arguably 



The experimental limits on the decay rates of the exclusive decays B 
being searched for by the CLEO and CDF collaborations [0, 110 , 144 
closer to the SM-based estimates, can only be interpreted in specific models of form 
factors, which hinders somewhat their transcription in terms of the information on the 
underlying Wilson coefficients. However, many of these form factors can be related to the 
known ones using ideas from heavy quark effective theory augmented with QCD sum rule 
estimates in studying the heavy to light form factors in B decays [ 145 |. We will not take 
up such exclusive decays, important as they are, in these lectures and content ourselves 
with presenting the expected branching ratios for some of the experimentally interesting 
decay modes. Using the exclusive-to-inclusive ratios Rru = ^{B Ki'^i")/T{B 
XJ+£-) = 0.07 ± 0.02 and = T{B K*i+i-)/T{B XJ+i') = 0.27 ± 0.0.07, 

which were estimated in[ 102|] , the results are presented in Table 0. 

In conclusion, the semileptonic FCNC decays B — >■ Xsi~^i~ (and also the exclusive 
decays) will provide very precise tests of the SM, as they will determine the signs and 
magnitudes of the three Wilson coefficients, C7, Cg-^ and Ciq. This, perhaps, may also 
reveal physics beyond-the-SM if it is associated with not too high a scale. The MSSM 
model is a good case study where measurable deviations from the SM are anticipated and 



worked out 129 



3.9 Summary and overview of rare B decays in the SM 

The rare B decay mode B — ^ X^i^z/, and some of the exclusive channels associated with 
it, have comparatively larger branching ratios. The estimated inclusive branching ratio 



58 



in the SM is M, 



B{B Xsvv) = (4.0 ± 1.0) X 10" 



(218) 



where the main uncertainty in the rates is due to the top quark mass, as the top quark 
contribution completely dominates the decay rate, and a residual one from the semilep- 
tonic branching ratio Bsl- The scale-dependence, which enters indirectly through the top 
quark mass, has been brought under control through the NLL corrections, calculated in 

Xdi^u is related by the ratio of the 



||147|| . The corresponding CKM-suppressed decay B 
CKM matrix element squared [|102|| : 



B{B Xsuu) 



td\ 



Wt 



ts 



(219) 



Similar relations hold for the ratios of the exclusive decay rates, in which the r.h.s. 
depends additionally on the ratios of the form factors squared, which deviate from unity 
through S'f/(3)-breaking terms, in close analogy with the exclusive radiative decays dis- 
cussed earlier. These decays are particularly attractive probes of the short- distance 
physics, as the long-distance contributions are practically absent in such decays. Hence, 
relations such as the one in ( |219D provide, in principle, one of the best methods for the 
determination of the CKM matrix element ratio [|102|| . From the practical point 

of view, however, these decay modes are rather difficult to measure, in particular at the 
hadron colliders and probably also at the B factories. The best chances are at the 
decays at LEP, and indeed the present best upper limit is based on the analysis of LEP 
data ||146|| . This derived limit ||146| 

(220) 



B{B Xvu) < 3.9 X 10" 



will be hard to improve significantly in the foreseeable future. The estimated branching 
ratios in a number of inclusive and exclusive decay modes are given in Table ^ updating 
the estimates in [ |102|] . 

Further down the entries in Table ^ are listed some two-body rare decays, such as 
studied in in ||150|| - |p.54|| , where only the lowest order contributions are 

?°,S0) ^ studied 

Some of them, 



{BIB',) 



77, 



calculated, i.e., without any QCD corrections, and the decays {B 
in the next-to-leading order QCD in 



147 



B^ — > /i+zi and perhaps also the radiative decay S° 



in particular, the decays 
77, have a fighting chance 



to be measured at LHC. The estimated decay rates, which depend on the pseudoscalar 
coupling constant /^^ (for ^B^-decays) and fs^ (for .B^-decays), together with the present 
experimental bounds are listed in Table |^. Since no QCD corrections have been included 
in the rate estimates of {Bg, Bet) 11, the branching ratios are somewhat uncertain. The 
constraints on beyond-the-SM physics that will eventually follow from these decays are 
qualitatively similar to the ones that (would) follow from the decays B ^ X^ + '~f and 
B Xsi~^i~, which we have discussed at length earlier. 



4 An Update of the CKM Matrix 



In updating the CKM matrix elements, the Wolfenstein parametrization [jl4| has been 
used, which has been given earlier. The emphasis here is on those quantities which 
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constrain the CKM parameters, A, A, p and rj. However, for the sake of completeness we 
also quote the values for the other elements of the CKM matrix. 

We recall that \ Vus \ has been extracted with good accuracy from K irev and hyperon 
decays to be 

|V;,| = A = 0.2205 ±0.0018 . (221) 
This agrees quite well with the determination of Vud — 1 ~ ^-^^ from /3-decay, 



IK 



ud\ 



0.9744 ±0.0010 . 



The values of the matrix elements involving the charm quark row are: 



1.01 ±0.18 ; 



cd 



0.204 ±0.017 



(222) 



(223) 



The value of the third matrix element involving the charm quark Vcb is rather important 
for the CKM fits, as it determines the parameter A. This has been discussed earlier, 
yielding: 

iKfel = 0.0388 ±0.0036 ^ A = 0.80 ± 0.075 . (224) 

The other two CKM parameters p and rj are constrained by the measurements of 
iKife/Kfel, |e| (the CP-violating parameter in the kaon system), Xd {B^-B^ mixing) and (in 
principle) e'/e (AS" = 1 CP- violation in the kaon system). The constraints from e'/e are 
not included, due to the various experimental and theoretical uncertainties surrounding 
it at present. For an up to date review of this topic, we refer to [^, which contains an 
exhaustive list of references to the original literature on this and other related topics. 

The ratio iKih/Kbl is determined to be: 



V,. 



ub 



cb 



0.08 ±0.02 



p2 + ^2 = 0.36 ±0.08 . 



The experimental value of |e| is 



(2.26 ±0.02) X 10" 



(225) 



(226) 



Theoretically, |e| is essentially proportional to the imaginary part of the box diagram for 
K^-K'^ mixing and is given by ||155|| 



GlflMKMl 



Bk [A^X^rj] (yc {VcifsiVc, Vt) - Vcc] 



+ Vuytf2{yt)A'X' 



P)), 



(227) 



where i/i = m'f/M^, and the functions /2 and /s can be found in Ref. ||120|| . Here, the fj. 



are QCD correction factors, of which fjcc ||156|| and fju [ P-57|| were calculated some time ago 
to next-to- leading order, and fjct was known only to leading order ||158| , |159|| . Recently, 
this last renormalization constant was also calculated to next-to- leading order [|160|| . In 



the CKM fits published in [83|, the following values for the renormalization-scale- invariant 



coefficients have been used: fjcc — 1.32, rju — 0.57, fjct — 0.47, calculated for rhc = 1.3 



GeV and the NLO QCD parameter A^^^ = 310 MeV in Ref. [|T6 
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The final parameter in the expression for |e| is the renormalization-scale independent 
parameter Bk, which represents our ignorance of the hadronic matrix element 
{K^\{d'y'^{l — 75)5) \K^). The evaluation of this matrix element has been the subject of 
much work. The earlier results are summarized in Ref . [[L61|| . Some recent calculations of 
I3x using the lattice QCD methods and 1/Nc approach ||165 



are: Bk = 0.83 ± 0.03 
= 0.67 ±0.07 [JLQCD 
and Bk = 0.70 ±0.10 



[Sharpe ||T62[], Bk = 0.86 ± 0.15 [APE Collaboration ^6§], Bk 
Collaboration [|l^], Bk = 0.78 ± 0.11 [Bernard and Soni |164 
[Bijnens and Prades [p.65||]. 



We now turn to B^-B^ mixing. 



The present world average of = AM^/Frf, which is 



a measure of this mixing, is |167 



Xd = 0.71 ±0.04 



(228) 



which is based on time-integrated measurements which directly measure Xd, and on time- 
dependent measurements which measure the mass difference AMd directly. This is then 
converted to Xd using the i?^ lifetime, which is known very precisely [r^Bd) = 1.56 ± 
0.05 ps) [||]. 

From a theoretical point of view it is better to use the mass difference AMd, as it 
liberates one from the errors on the lifetime measurement. In fact, the present precision 
on AMd, pioneered by time-dependent techniques at LEP, is quite competitive with the 
precision on Xd- The LEP average for AMd has been combined with the one from CDF 
and that derived from time-integrated measurements yielding the present world average 
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AMd = 0.457 ±0.019 (ps)" 



(229) 
Unlike the 

kaon system, where the contributions of both the c- and the t-quarks in the loop were 
important, this diagram is dominated by t-quark exchange: 



The mass difference AMd is calculated from the B^-B^ box diagram. 



AMd 



{fBBB,)fjBytf2{yt)\v:dVtt\ 



(230) 



where, using Eq. 
In Ref, 



KdVttl 



Here, 17^ is the QCD correction. 



157|| , this correction is analyzed including the effects of a heavy t-quark. It is 



found that fis depends sensitively on the definition of the t-quark mass, and that, strictly 
speaking, only the product fjB{yt)f2{yt) is free of this dependence. In the fits presented 
here we use the value fjB = 0.55, calculated in the MS scheme, following Ref. ||157| . 
Consistency requires that the top quark mass be rescaled from its pole (mass) value to 



the value mj(mi(pole)) in the MS scheme, which is typically about 10 GeV smaller |6^. 



For the B system, the hadronic uncertainty is given by Bb^, analogous to Bk in 



the kaon system, except that in this case, also /^^ is not measured. In the CKM fits [5^ 



the following ranges for Bb^ and Bb^, which are compatible with results from both 
lattice-QCD m, W 



and QCD sum rules ||168|| , were used 



180 ± 50 MeV 



Bb, = 1.0 ±0.2. 



(231) 
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With this theoretical input, one can get an estimate of the CKM matrix element \Vtd\- 
With AMd known very precisely, and Arrit/mt = ±6% p|, the error on \ Vtd\ is dominated 
by theoretical error: 

\Vtd\ = (0.92 ± 0.02 ± 0.10 ± 0.28) x 10"^ . (232) 

The three errors are from AMd (expt), Arrit and theory, respectively. Adding the errors 
in quadrature, this gives 

\Vtd\ = (0.92 ± 0.30) X 10~2 , (233) 

which is better than the range 0.004 < \Vtd\ < 0.015 following from unitarity 0]. 

To complete the estimates of all the nine CKM matrix elements, we list here the 
determination of | VJs|/| Kbl from the inclusive branching ratio B{B — >■ + 7), given in 

0.85±0.22(expt + th), (234) 



Eq. (M: 

\Vts\ 



\Vcb\ 

where, like in \ Vcb\, we have added the experimental and theoretical errors linearly. Com- 
bining it with \V^b\ = 0.0388 ± 0.0036, gives 

\Vts\ = 0.033 ± 0.009(expt + th). (235) 



A determination of \Vtb\ can be derived from the CDF measurements Eq. (p!9|). While 
compatible with the unitarity bound 

0.9988 < \Vtb\ < 0.9995 , (236) 

the direct determination of \Vtb\ from CDF is far less accurate. The present knowledge of 
the nine CKM matrix elements is summarized in Table |^. 

4.1 The present profile of the Unitarity Triangle 

The entries in Table ^ provide a test of unitarity. Since, the value quoted for \Vtb\ is from 
unitarity, the other two constraints are: 

iKdl' + lK^ + lKbP = 0.998 ±0.002, 

\Vcd\^ + \Vcs\^ + \Vcb\^ = 1.03 ±0.18. (237) 

This shows that unitarity is well satisfied. One could now ask the question, how well are 
the parameters p and f] determined at present, which in turn determine the profile of the 
unitarity triangles. This is done by fitting the CKM parameters using the quantities we 
have discussed. 

In order to find the allowed unitarity triangles, the computer program MINUIT was 
used to fit the CKM parameters A, p and rj to the experimental values of \Vcb\, \Vub/Vcb\, 
\e\ and Xd- Since A is very well measured, it was fixed to its central value given above. 



Two types of fits can be attempted |8^ 



Fit 1: the "experimental fit." Here, only the experimentally measured numbers are 



used as inputs to the fit with Gaussian errors; the coupling constants fB^yBs^ and 
Bk are given fixed values. 
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• Fit 2: the "combined fit." Here, both the experimental and theoretical numbers are 
used as inputs assuming Gaussian errors for the theoretical quantities. 

We first discuss the "experimental fit" (Fit 1). The goal here is to restrict the allowed 



range of the parameters (p, r]) for given values of the coupling constants /b^v-Bb^ and 



Bk- For each value of and fB^yBB^, the CKM parameters A, p and rj are fit to the 
experimental numbers given earlier and the is calculated. In the CKM fits performed 
in ||8^, specific values in the range 0.4 to 1.0 for B^ were considered and it was shown 
that for Bk = 0.4 a very poor fit to the data is obtained, so that such small values are 
quite disfavoured. The fits performed in the range 0.6 < Bk < 1.0, which adequately 
cover the more recent predictions given above, had a good quality. 

In Fit 2 in ||8^, a central value plus an error to Bk was assigned and two ranges for 
Bx were considered: 

Bk = 0.8 ± 0.2 , (238) 
which reflects the estimates of this quantity in lattice QCD, or 

Bk = 0.6 ± 0.2 , (239) 

which overlaps with the values suggested by the earlier chiral perturbation theory esti- 



mates 166 . It was shown that there is not an enormous difference in the results for the 



two ranges. However, as now there seems to be a theoretical consensus emerging, with 
Bk — 0.8, we show the case where Jb^ \J Bb^ is varied in the range 130 MeV to 230 MeV. 
The fits are presented as an allowed region in p-77 space at 95% C.L. (x^ = Xmm + 6-0). 
The results are shown in Fig. |ll]. As we pass from Fig. |ll](a) to Fig. |ll|(e), the unitarity 
triangles represented by these graphs become more and more obtuse. Even more striking 
than this, however, is the fact that the range of possibilities for these triangles is quite 
large. There are two things to be learned from this. First, our knowledge of the unitarity 
triangle is at present rather poor. Second, unless our knowledge of hadronic matrix ele- 
ments improves considerably, measurements of |e| and Xd-, no matter how precise, will not 
help much in further constraining the unitarity triangle. This is why measurements of 
CP- violating rate asymmetries in the B system are so important ||169| , [16| . Being largely 
independent of theoretical uncertainties, they will allow us to accurately pin down the 
unitarity triangle. With this knowledge, we could deduce the correct values of ]3k and 



Jb^^Bb^, and thus rule out or confirm different theoretical approaches to calculating 
these hadronic quantities. 



Despite the large allowed region in the p-r] plane, certain values of Bk and fB^\j Bb^ 
are disfavoured since they do not provide a good fit to the data. For example, fixing 
I3k = 0.8, which is presently theoretically favoured, we can use the fitting program to 



provide the minimum for various values of fB^y Bb^- The results are shown in Table 
along with the best fit values of {p,Ti). Since we have two variables (p and 77), we 



use Xmin < 2.0 as our "good fit" criterion, and we see that /b^v-Bb^ < 130 MeV and 



fBj^y Bbj > 240 MeV give poor fits to the existing data. 

We now discuss the "combined fit" (Fit 2). Since the coupling constants are not known 
and the best we have are estimates given in the range in Eq. ( p38| ), a reasonable profile of 
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Figure 11: Allowed region in p-rj space, from a fit to the experimental values given in 
the text. We have fixed Bk = 0.8 and vary the coupling constant product /^^ \J Bb^ 
as indicated on the figures. The solid line represents the region with = Xmin + 6 
corresponding to the 95% CL. region. The triangles show the best fit. (Figure taken 
from H.) 
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Figure 12: Allowed region in p-rj space, from a simultaneous fit to both the experimental 
and theoretical quantities discussed in the text. The theoretical errors are treated as 
Gaussian for this fit. The solid line represents the region with = Xmin~^^ corresponding 
to the 95% CL. region. The triangle shows the best fit. (Figure taken from |83|.) 



the unitarity triangle at present can be obtained by letting the coupling constants vary in 
this range. The resulting CKM triangle region is shown in Fig. As is clear from this 
figure, the allowed region is rather large at present. The preferred values obtained from 
the "combined fit" are 

(p, n) = (-0.07, 0.34) (with = 6.6 x 10^^) . (240) 



4.2 Xs and the Unitarity Triangle 



Mixing in the B^-B'^ system is quite similar to that in the B^-B^ system. The B^^-B^ box 
diagram is again dominated by t-quark exchange, and the mass difference between the 
mass eigenstates AM^ is given by a formula analogous to that of Eq. 



(^2 



(241) 



Using the fact that \Vch\ = \Vts\, it is clear that one of the sides of the unitarity triangle, 
|V^d/AV^{,|, can be obtained from the ratio of AM^^ and AM^, 



AM, fiBMBAfiBB^ 



AM, f^^^MB, {PbBb^ 



ts 



V, 



td 



(242) 



All dependence on the t-quark mass drops out, leaving the square of the ratio of CKM 
matrix elements, multiplied by a factor which reflects 5"^/ (3) flavour breaking effects. The 
only real uncertainty in this factor is the ratio of hadronic matrix elements. Whether 
or not Xs can be used to help constrain the unitarity triangle will depend crucially on 
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the theoretical status of the ratio Jb^BbJ/b Bb^- In what follows, we will take = 

UBs4^s)l{fBa4^d) = (1-16 ±0.1), consistent with both lattice-QCD [|4[ |32| and 
QCD sum rules ^6^. (The SU(3)-breaking factor in AM./AM^ is ^l) 

The mass and lifetime of the Bg meson have now been measured at LEP and Tevatron 
and their present values are Mb^ = 5370.0 ± 2.0 MeV and t{Bs) = 1.55 ± 0.10 ps ||69| . 



The QCD correction factor fjB^ is equal to its counterpart, i.e. fjB^ = 0.55. The main 
uncertainty in Xg (or, equivalently, AMg) is now Jb^Bb,- Using the determination of A 
given previously, tb, = 1.55 ± 0.10 (ps) and mj = 170 ±11 GeV, we obtain 

AM, = (13.1 ± 2.8) —^^^^^ (ps)-\ 
^ ' (230 MeV)2 ' ' 

Xs = (20.3 ± 4.5) — ^^^—^ . (243) 
^ ' (230 MeV)2 ^ ' 



The choice /bs^Bbs = 230 MeV corresponds to the central value given by the lattice- 
QCD estimates, and with this the CKM-fits discussed earlier give Xg — 20 as the preferred 
value in the SM. Allowing the coefficient to vary by ±2a, and taking the central value for 



fBsyBBs, this gives 

11.4 < Xs< 29.4 , 
7.5 (ps)-^ < AM, < 18.7 {ps)~^ . (244) 

It is difficult to ascribe a confidence level to this range due to the dependence on the un- 



known coupling constant factor. In particular, the upper limit is scaling as [/^^ y i?Bs/(230MeV)]^ 
All one can say is that the standard model predicts large values for Xg, most of which are 
above the present experimental limit Xg > 8.8 (equivalently AM, > 6.1 {ps)~^) ||126| , |167 . 



The ALEPH lower bound AM, > 6.1 (ps)"^ (95% C.L.) [g and the present world 
average AM^ = (0.457±0.019) (ps)^^ can be used to put a bound on the ratio AM^/AM^. 
The lower limit on AM, is correlated with the value of /,, the fraction of b quark frag- 
menting into Bg meson, as shown in the ALEPH analysis [|126|| . The value obtained 
from the measurement of the quantity fgBR{Bg Dgiui) is /, = (11.0 ± 2.8)%. The 
time-integrated mixing ratios x ^"^^ Xd, assuming maximal mixing in the Bg-Bg sys- 
tem Xs = 0.5, give /, = (9.9 ± 1.9)%. The weighted average of these numbers is 
/, = (10.2±1.6)% [le^. With /, = 10%, one gets AM, > 5.6 {ps)~^ at 95% C.L., yielding 
AMg/AMd > 11.8 at 95% C.L. Assuming, however, /, = 12% gives AM, > 6.1 {ps)-\ 
yielding AMg/AMd > 12.8 at 95% C.L. We will use this latter number. 

The 95% confidence limit on AMg/AMd can be turned into a bound on the CKM 
parameter space (p, r]) by choosing a value for the SU(3)-breaking parameter We 
assume three representative values: = 1.1, 1.35 and 1.6, and display the resulting 
constraints in Fig. |T^. From this graph we see that the ALEPH bound marginally restricts 
the allowed p-rj region for small values of but does not provide any useful bounds for 
larger values. 

Summarizing the discussion on x,, we note that the lattice-QCD- inspired estimate 
fBs^BBs — 230 MeV and the CKM fit predict that Xg lies between 12 and 30, with a 
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Figure 13: Further constraints in p-r] space from the ALEPH bound on AM^. The bounds 
are presented for 3 choices of the SU(3)-breaking parameter: = 1.1 (dotted hne), 1.35 
(dashed hne) and 1.6 (sohd hne). In ah cases, the region to the left of the curve is ruled 
out. (Figure taken from [^.) 



central value around 20. The upper and lower bounds and the central value scale as 



(/ss V -^Bs/230 MeV)^. The present constraints from the lower bound on Xg on the CKM 
parameters are marginal but this would change with improved data. In particular, one 
expects to reach a sensitivity Xg ~ 15 (or AMg ~ 10 ps~^) at LEP combining all data 



and tagging techniques ||167|| . One expects comparable sensitivity at the SLC, where the 
beam polarization can be used advantageously, and at the HERA-B, CDF and DO where 
the Bs mesons will be produced copiously. The entire range for Xg for the SM given in 
Eq. ( p44| ) will be accessible at the LHC. A measurement of Xs (equivalently AM^) would 
be very helpful in further constraining the CKM parameter space. 



5 CP Violation in the B System 



This topic has been reviewed in [^, and since only marginal changes have taken place 
in the experimental situation and our present understanding of CP violation in the B 
system, we reproduce essentially this discussion here. 

It is expected that the B system will exhibit large CP-violating effects, characterized 
by nonzero values of the angles a, (3 and 7 in the unitarity triangle (Fig. [1|) [|169|| . The 
most promising method to measure CP violation is to look for an asymmetry between 
r(S° — s> /) and T{B^ /), where / is a CP eigenstate. If only one weak amplitude 
contributes to the decay, the CKM phases can be extracted cleanly (i.e. with no hadronic 

uncertainties). Thus, sin 2a, sin 2(3 and sin 27 can in principle be measured in B^ —>■ tt+tt", 

'-Brf J/ijjKs and^Bg pKs, respectively. 

Unfortunately, the situation is not that simple. In all of the above cases, in addition 
to the tree contribution, there is an additional amplitude due to penguin diagrams ||170|| . 
In general, this will introduce some hadronic uncertainty into an otherwise clean measure- 
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ment of the CKM phases. In the case of^Ba J/ipKs, the penguins do not cause any 
problems, since the weak phase of the penguin is the same as that of the tree contribution. 
Thus, the CP asymmetry in this decay still measures sin 2(3. 

For Bd — > TT+TT^, however, although the penguin is expected to be small with respect 
to the tree diagram, it will still introduce a theoretical uncertainty into the extraction 
of a. Fortunately, this uncertainty can be removed by the use of isospin | |1 71| | . The 
key observation is that the 1 = 2 component of the B ^ mr amplitude is pure tree 
(i.e., it has no penguin contribution) and therefore has a well-defined CKM phase. By 
measuring the rates for B~^ 7r"^7r°, B^ — >■ tt^tt" and B^ — >■ 7r°7r°, as well as their CP- 
conjugate counterparts, it is possible to isolate the 1 = 2 component and obtain a with 
no theoretical uncertainty. Thus, even in the presence of penguin diagrams, sin 2a can 
in principle be extracted from the decays B ^ tttt. It must be admitted, however, that 
this isospin program is ambitious experimentally. If it cannot be carried out, the error 
induced on sin 2a is of order \P/T\, where P (T) represents the penguin (tree) diagram. 
The ratio \P/T\ is difficult to estimate - it is dominated by hadronic physics. However, one 
ingredient is the ratio of the CKM elements of the two contributions: \V^lVtd/V*,jVud\ — 
\Vtd/Vub\- From the fits in |]S3|, the allowed range for the ratio of these CKM matrix 



elements is 

1.2 < 



ub 



< 5.8 , (245) 



with the central value close to 3. 

It is Bs —>■ pKs which is most affected by penguins. In fact, the penguin contribution 
is probably larger in this process than the tree contribution. This decay is clearly not 
dominated by one weak (tree) amplitude, and thus cannot be used as a clean probe of 
the angle 7. Instead, two other methods have been devised, not involving CP-eigenstate 
final states. The CP asymmetry in the decay Bg DfK^ can be used to extract sin^ 7 
172 |. Similarly, the CP asymmetry in B"^ D^pK^ also measures sin^7 [173|. Here, 



D^p is a or which is identified in a CP-eigenstate mode (e.g. tt+tt", K^K^ , ...). 

These CP-violating asymmetries can be expressed straightforwardly in terms of the 
CKM parameters p and 77. The 95% C.L. constraints on p and rj found previously can 
be used to predict the ranges of sin 2a, sin 2/5 and sin^7 allowed in the standard model. 
The allowed ranges which correspond to each of the figures in Fig. 11, obtained from Fit 

1, are found in Table Bl In this table we have assumed that the angle (3 is measured in 

( — ) 

Bd J/'^Ks, and have therefore included the extra minus sign due to the CP of the 
final state. 

Since the CP asymmetries all depend on p and rj, the ranges for sin 2a, sin 2/5 and 
sin^ 7 shown in Table ^ are correlated. That is, not all values in the ranges are allowed 
simultaneously. We illustrate this in Fig. |14|, corresponding to the "experimental fit" (Fit 
1), by showing the region in sin2a-sin2/? space allowed by the data, for various values 



of fBayBsa- Given a value for f^^yBB^, the CP asymmetries are fairly constrained. 



The parameters used in the central figure in Fig. |14|, namely /b^y-B_b^ = 180 MeV 
and Bk = 0.8, are the best theoretical estimates at present. This figure and the third 
row in Table ^ then yield: sin 2/3 > 0.42 and sin^7 > 0.68. However, since there is still 
considerable uncertainty in the values of the coupling constants, a more reliable profile of 
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the CP asymmetries at present is given by the "combined fit" (Fit 2) of [Q, where present 
theoretical and experimental values have been convoluted in their allowed ranges. The 
resulting correlation is shown in Fig. |15[ From this figure one sees that the smallest value 
of sin 2/? occurs in a small region of parameter space around sin 2a ~ 0.4-0.6. Excluding 

this small tail, one expects the CP-asymmetry in J/'^Ks to be at least 30%. 

It may be difficult to extract 7 using the techniques described above. First, since 
Bg DfK^ involves the decay of Bg mesons, such measurements must be done at 
hadron colliders. At present, it is still debatable whether this will be possible. Second, the 
method of using B^ — > D'^pK^ to obtain 7 requires measuring the rate for B~^ —>■ D^K^ . 
This latter process has an expected branching ratio of < O(10~^), so this too will be hard. 

Recently, a new method to measure 7 was proposed ||174|| . Using a flavour SU(3) 
symmetry, along with the neglect of exchange- and annihilation-type diagrams, it was 
suggested that 7 could be found by measuring rates for the decays B^ — > t[^K^ , B^ 
7r+i^°, B~^ tt+tt", and their charge-conjugate processes. The irK final states have both 
1=1/2 and 1 = 3/2 components. The crucial ingredient is that the gluon- mediated 
penguin diagram contributes only to the J = 1/2 final state. Thus, a linear combination 
of the B^ tt'^K'^ and B^ n^K^ amplitudes, corresponding to / = 3/2 in the irK 
system, can be related via flavour SU(3) to the purely 1 = 2 amplitude in B^ tt+tt", 
permitting the construction of an amplitude triangle. The difference in the phase of the 
B~^ — > 7r"'"7r° side and that of the corresponding triangle for B~ decays was found to be 
27. SU(3) breaking can be taken into account by including a factor fx/fn in relating 
i? — > TTTT decays to the B 71K decays ||175 . 



The key assumption is that the penguin is mediated by gluon exchange. However, 
there are also electroweak contributions to the penguins [|176|| . These electroweak penguins 
(EWP's) are not constrained to be isosinglets. Thus, in the presence of EWP's, there is 
no longer a triangle relation B ttK and B ^ tttt amplitudes [|177|| . Indeed, electroweak 
penguins can, in principle, even invalidate the isospin analysis in B ^ tttt, since the 
1 = 2 amplitude will include a contribution from EWP's, and hence will no longer have 
a well-defined weak CKM phase. However, theoretical estimates [ p,77| , p,78|| show that 
electroweak penguins are expected to be relatively unimportant for B ^ -k-k. 

The question of the size of EWP's has therefore become a rather interesting question, 
and a number of papers have recently appeared discussing this issue ||179|| . These include 
both theoretical predictions, as well as ways of isolating EWP's experimentally. The 
general consensus is that EWP's are large enough to invalidate the method of Ref. ||174|| for 
obtaining 7. However, two new methods making use of the flavour SU(3) symmetry, and 
which do not have any problems with electroweak penguins, have been suggested. Both 
are rather complicated, making use of the isospin quadrangle relation among B ttK 
decays, as well as B~^ tt+tt" plus an additional decay: Bs rin^ in one case ||178|| , 
5+ —>■ TjK^ with 7] = 7]^ m the other [|180|| . Although electroweak penguins do not 
cause problems, SU(3)-breaking effects which cannot be parametrized simply as a ratio 
of decay constants are likely to introduce errors of about 25% into both methods. It 
is clear that this is a subject of great interest at the moment, and work will no doubt 
continue. Recently, another quadrangle relation involving charged B decays B^ —>■ rjK^ 
and B'^ — > r]'K^ along with the decays B~^ — > n^K^ and B^ 7r^K~^ has been proposed 
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Figure 14: Allowed region of the CP asymmetries sin 2a and sin 2/5 resulting from the 
"experimental fit" of the data for different values of the coupling constant fs^ \/^b^ 
indicated on the figures a) - e). We fix Bk = 0.8. (Figure taken from [Q.) 
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= 180 + 50 MeV , B„ = 1 .0 + 0.2 



= 0.8 + 0.2 



-0.4 -0.2 0.2 
sin 2a 



Figure 15: Allowed region of the CP asymmetries sin 2a and sin 2/3 resulting from the 



"combined fit" of the data for the ranges for fB^y Bb^ ^-nd Bk given in the text. (Figure 
taken from 



to determine the weak phase 7 

5.1 Summary of the CKM fits and CP asymmetries in B Decays 

We summarize the results of this section: 

(i) We have presented an update of the CKM unitarity triangle using the theoretical 
and experimental improvements in the following quantities: \Vcb\, \Vub/Vcb\, AM^, T{Bd), 
Vcc, Vet- The fits can be used to exclude extreme values of the pseudoscalar coupling 



constants, with the range 130 MeV < fB^yBBa < 270 MeV still allowed for Bk = 1. The 
lower limit of this range is quite i?i^-independent, but the upper limit is strongly correlated 



with the value chosen for Bk- For example, for Bk = 0.8 and 0.6, fB^yBB^ < 240 and 
210 MeV, respectively, is required for a good fit. The solutions for ]3k = 0.8 ± 0.2 are 
slightly favoured by the data as compared to the lower values. These numbers are in very 
comfortable agreement with QCD-based estimates from sum rules and lattice techniques, 
and recently also with estimates from improved chiral perturbation theory. The statistical 
significance of the fit is, however, not good enough to determine the coupling constant 



more precisely. The fits in ||8^ show that Bk < 0.4 is strongly disfavoured by the data, 
since the quality of fit for such values is very poor. This value of ]3k is now of very little 
theoretical interest, anyway. 

(ii) The newest experimental and theoretical numbers restrict the allowed CKM uni- 
tarity triangle in the (p, ?7)-space somewhat more than before. However, the present 
uncertainties are still enormous - despite the new, more accurate experimental data, our 
knowledge of the unitarity triangle is still poor. This underscores the importance of mea- 
suring CP-violating rate asymmetries in the B system. Such asymmetries are largely 
independent of theoretical hadronic uncertainties, so that their measurement will allow 
us to accurately pin down the parameters of the CKM matrix. Furthermore, unless our 
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knowledge of the pseudoscalar coupling constants improves considerably, better measure- 
ments of such quantities as Xd will not help much in constraining the unitarity triangle. 
On this point, help may come from the experimental front. It may be possible to measure 
the parameter Z^^, using isospin symmetry, via the charged- current decay t^u^-. 
With iKfe/Kfel = 0.08 ± 0.02 and /b, = 180 ± 50 MeV, one gets a branching ratio 
BR{B^ T^u^) = (1.5-14.0) X 10-^ with a central value of 5.2 x 10"^ This lies in the 
range of the future LEP and asymmetric 5-factory experiments, though at LEP the rate 
Z — i> BcX — >■ T^UrX could be just as large as Z ^ B^X r^UrX. Along the same 
lines, the prospects for measuring {fBj_,fBs) in the FCNC leptonic and photonic decays 
of B^ and 5° hadrons, {B^, 5°) /U^/i", (-8°, B^) 77 in future B physics facilities are 
not entirely dismal. 

(iii) We have determined bounds on the ratio |Vtd/Vis| from our fits. For 130 MeV < 
Jb^ \J Bb^ < 270 MeV, i.e. in the entire allowed domain, at 95 % CL. we find 

Vtd 



0.13 < 



ts 



< 0.35 



(246) 



The upper bound from our analysis is more restrictive than the current experimental upper 
limit following from the CKM-suppressed radiative penguin decays BR{B ^ + 7) and 



BR{B ^ p + -f), which at present yield at 90% C.L. |4 

Vtd 



Vt 



ts 



< 0.64-0.75 



(247) 



depending on the model used for the SU(3)-breaking in the relevant form factors \\121[ |122 . 
Long-distance effects in the decay B^ + 7 may introduce theoretical uncertainties 

comparable to those in the SU(3)-breaking part but the corresponding effects in the decays 



B^ {p^,uj) + 7 are expected to be very small ||124|| . Furthermore, the upper bound is 



now as good as that obtained from unitarity, which gives 0.08 < |Vtd/Vts| < 0.36, but the 
lower bound from our fit is more restrictive, 
(iv) Using the measured value of uit, we find 



(20.3 ±4.5) 



(230 MeV)^ 



(248) 



Taking fBsyBBs = 230 (the central value of lattice-QCD estimates), and allowing the 
coefficient to vary by ±2a, this gives 



11.4 < X, < 29.4 . 



(249) 



No reliable confidence level can be assigned to this range - all that one can conclude is 
that the SM predicts large values for Xg, most of which lie above the ALEPH 95% C.L. 
lower limit of Xg > 8.8. 

(v) The ranges for the CP- violating rate asymmetries parametrized by sin2Q;, sin 2/3 
and and sin^7 are determined at 95% C.L. to be 



-1.0 < sin 2a < 1.0 , 
0.21 < sin 2/? < 0.93 , 
0.12 < sin^-f < 1.0 . 



(250) 
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(For sin 2a < 0.4, we find sin 2/3 > 0.3.) Electroweak penguins may play a significant role 
in some methods of extracting 7. Their magnitude, relative to the tree contribution, is 
therefore of some importance. One factor in determining this relative size is the ratio of 
CKM matrix elements |Vtrf/Kife|. We find 

< 5.8 . (251) 



1.2 < 



Vt 



td 



n.h 
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f R , R ^ ^ X T+r" 


f2 fi + n 5") X If)-'' 




f R , R ^ ^ X/T+r" 


fl 5+1 X If)-^ 

I X . tJ _l X . / /N X w 




f R , R ^ ^ Ke+e^ 


('5 9 + 2 3) X 10"'' 


< 1.2 X 10-5 CLEO g 


( n, R ^ _^ Kii+ir 
yj-'di J-'u) ' i\ i-i ft 


n + 1 5'i X 1 n^'^ 


< 0.9 X 10-5 CLEO H 


fR. R ^ K*P^P~ 


('2 3 + n Q"! X 1 


< 1.6 X 10-5 CLEO g 


fR, R ) ^ K*U + U~ 
y^di ^u) ' fi fi 


fi 5 + n fi^i X in^^ 


< 2.5 X 10-5 CDF |148| 


f 5 , R ) X. Z/Z/ 

y^d-i-^uj ^^s 


(4.0 ± 1.0) X 10"^ 


< 3.9 X 10-^ |ig 


y^di ^u) ' ^d '^'^ 


f2 3 + 2 0^1 X 10"^ 




y-'-'di -'-'uj ' -f^ 1^1^ 


f3 2 + 1 fi") X 10~^ 






n 1 + 55"! X 1 0^^ 




J-Js J J 


f3 ± 1 0") X 10"'' 


< 1.1 X 10-^ L3 |: 


L4 


9| 


Bd 77 


(1.2 + 1.1) X 10-^ 


< 3.8 X 10-5 L3 |149| 


fi, ^ r+r- 


(7.4 + 2.1) X 10-^ 




fid ^ r+r- 


(3.1 + 2.9) X 10-^ 




Bs fi^fi' 


(3.5 + 1.0) X 10-9 


< 8.4 X 10-^ CDF 


148| 


Bd 


(1.5 + 1.4) X 10-1° 


< 1.6 X 10-^ CDF 


148| 


Bs e+e- 


(8.0 + 3.5) X 10-1^ 




5d ^ e+e- 


(3.4 + 3.1) X 10-1^ 





Table 2: Estimates of the branching fractions for FCNC -B-decays in the standard model 
taking into account the uncertainties in the input parameters as discussed in the text. The 
entries in the second column correspond to the short-distance contributions only except 
for the radiative decays 5„ — + 7 and Bd — >■ [p^^oj) + 7, where long-distance effects 
have also been included. For the two-body branching ratios, we have used /b^ = 200 MeV 
and Jbs/ = 1-16. The CKM matrix element ratio is taken as |Vf(i|/|Vts| = 0.24 + 0.11, 
as determined from the present fits, and this error is folded in quadrature with the other 
errors in estimating the relevant branching fractions. Experimental measurements and 
upper limits are also listed. Note that the limit quoted from [|146|| is an indirect one. 
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1 T 7" 1 

\Vij\ 


rresent Value 




\Vud\ 


0.9744 ± 0.0010 [|] 




IK.I 


0.2205 ± 0.0011 [|] 




\Vub\ 


(3.1 ± 0.8) X 10 




\Vcd\ 


0.204 ±0.017 § 




\Vcs\ 


1.01 ±0.18 § 




\Vcb\ 


0.0388 ± 0.0036 [0 




\Vtd\ 


(9.2 ±3.0) X 10-3 




\Vts\ 


0.033 ±0.009 




\Vtb\ 


0.9991 ± 0.0004 [|] 





Table 3: Present values of the CKM matrix elements \Vij\. Note that the value of \Vtb\ 
follows from unitarity. All others are measured in decays discussed in the text and in the 
PDG review 0, from which references to the original literatue can be traced. 



/B,v/S^(MeV) 


{p,v) 


Anwn 


120 


(-0.42, 0.16) 


3.04 


130 


(-0.37, 0.19) 


1.32 


140 


(-0.32, 0.23) 


0.43 


150 


(-0.27, 0.26) 


0.07 


160 


(-0.22, 0.29) 


1.4 X 10-3 


170 


(-0.15, 0.32) 


0.05 


180 


(-0.08, 0.34) 


0.09 


190 


(-0.01, 0.35) 


0.06 


200 


(0.06, 0.35) 


0.01 


210 


(0.13, 0.35) 


0.02 


220 


(0.18, 0.35) 


0.2 


230 


(0.23, 0.35) 


0.61 


240 


(0.28, 0.35) 


1.29 


250 


(0.32, 0.35) 


2.22 



Table 4: The "best values" of the CKM parameters (p, rj) as a function of the coupling 
constant fs^ \J Bb^, obtained by a minimum fit to the experimental data, including the 
renormalized value of = 170 ± 11 GeV. We fix = 0.8. The resulting minimum 
values from the MINUIT fits are also given. (Table taken from [p3| .) 
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/B,^/S^(MeV) 


sin 2a 


sin 2/3 


sin^ 7 


130 
155 
180 
205 
230 


0.46 - 0.88 
0.75 - 1.0 
-0.59 - 1.0 
-0.96 - 0.92 
-0.98 - 0.6 


0.21 - 0.37 
0.31 - 0.56 
0.42 - 0.73 
0.49 - 0.86 
0.57-0.93 


0.12 - 0.39 
0.34 - 0.92 
0.68 - 1.0 
0.37- 1.0 
0.28 - 0.97 



Table 5: The allowed ranges for the CP asymmetries sin 2a, sin 2/5 and sin^7, corre- 
sponding to the constraints on p and 1] shown in Fig. Values of the coupling constant 

/b^ J are stated. We fix Bk = 0.8. The range for sin 2/5 includes an additional minus 



sign due to the CP of the final state J/^Ks- (Table taken from p3| 
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